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EXECUTIVE  SUMMARY 

In  1985,  the  Medical  Research  and  Evaluation  Facility  (MREF) 
developed  a  standardized  first-stage  screen  (MREF  Protocol  21,  May  1985)  to 
compare  liquid  or  powder  experimental  decontaminants  against  the 
dual -component  M258A1  skin  decontamination  kit  for  their  effectiveness  in 
mitigating  the  toxic  effects  of  percutaneous  exposure  to  organophosphate 
chemical  surety  materiels  (CSM).  The  testing  protocol  calls  for  strict 
standardization  of  methods,  materials,  and  agent  doses  across  the  individual 
screening  tests  that  are  spread  out  over  a  multi-year  time  frame.  A 
concurrent  standard  decontaminant  group  is  included  in  each  individual  test  of 
the  experimental  decontaminants.  Thus,  a  considerable  data  base  is  amassed 
over  time  pertaining  to  the  standard  decontamination  procedure  results. 

A  principal  objective  of  this  report  is  to  describe  and  illustrate 
statistical  methods  for  the  incorporation  of  the  historical  data  accumulated 
on  the  standard  decontaminant  results  to  enhance  the  statistical  sensitivity 
of  individual  comparisons  between  the  standard  and  experimental 
decontaminants.  The  basic  idea  is  that  the  historical  levels  and  variability 
of  the  standard  decontaminant  test  results  can  be  used  to  predict  a  likely 
range  for  the  concurrent  standard  decontaminant  test  results.  This 
information  can  be  incorporated  into  the  concurrent  test  procedures.  The 
tradeoff  in  using  this  historical  information  is  that  the  test  procedure  will 
be  more  sensitive  if  concurrent  results  fall  within  the  range  of  the  past 
results,  but  may  perform  worse  in  terms  of  the  Type  1  error  being  too  large  or 
too  small  than  a  test  that  ignores  the  historical  information  if  the 
.concurrent  standard  decontaminant  response  level  is  substantially  discrepant 
from  the  distribution  of  historical  response  levels.  For  this  reason,  the 
test  procedure  recommended  in  this  report  compares  the  experimental 
decontaminant  response  rate  to  a  weighted  average  of  the  concurrent  standard 
decontaminant  response  rate  and  the  historical  rate.  The  weight  associated 
with  the  historical  rate  increases  as  the  observed  time  to  time  variability  in 
the  historical  data  decreases  and  as  the  agreement  between  the  concurrent  rate 
and  the  historical  average  rate  increases. 


i  i 

Several  additional  statistical  aspects  of  the  screen  are  also 
discussed.  Control  chart  procedures  are  suggested  to  detect  drifts  over  time 
or  sudden  jumps  in  the  standard  decontaminant  responses  rates.  In  the  event 
the  control  charts  indicate  that  the  LD5I  for  the  standard  decontaminant  has 
shifted,  procedures  are  presented  for  carrying  out  studies  to  update  the  LDS0. 
These  procedures  are  designed  to  conserve  experimental  resources.  The  dose 
allocation  for  the  LDS(  studies  is  carried  out  in  a  stagewise,  adaptive 
fashion.  The  dose  selection  for  each  stage  of  the  design  is  based  on  the  test 
results  from  all  previous  stages.  It  is  designed  to  accommodate  unanticipated 
aspects  of  the  dose-response  relation.  The  stagewise,  adaptive  dose- 
allocation  strategy ' introduces  a  number  of  nonstandard  considerations  that 
necessitate  the  use  of  specialized  dose-response  model  fitting  procedures. 
Specialized  probit  analysis  model  fitting  procedures,  based  on  nonlinear 
regression  analysis,  are  discussed  and  illustrated  by  example. 
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LETHALITY  RATE  ESTIMATION  AND  TESTING  PROCEDURES 


1.0  INTRODUCTION 


In  1985  the  Medical  Research  and  Evaluation  Facility  (MREF) 
developed  a  standardized  first  stage  screening  test  entitled  "Assessment  of 
Liquid  or  Powder  Decontaminants  Against  GD,  Thickened  GD,  and  VX  Administered 
Topically  to  Rabbits"  (MREF  Protocol  21,  May  1985)  to  compare  liquid  or  powder 
experimental  decontaminants  against  the  dual -component  M258A1  skin 
decontamination  kit  for  their  effectiveness  in  mitigating  the  toxic  effects  of 
percutaneous  exposure  to  organophosphate  (OP)  chemical  surety  materiel'.:  (CSM) . 
The  standardized  screen  is  based  on  a  lethality  endpoint  in  laboratory  albino 
rabbits.  An  essential  aspect  of  this  testing  protocol  is  the  strict 
standardization  of  methods,  materials,  and  agent  dose  levels  (LDSI  doses 
associated  with  standard  decontamination  procedure)  that  are  used  to  screen 
numerous  experimental  decontaminants  throughout  a  time  period  extending  over 
multiple  years. 

Each  test  decontaminant  is  necessarily  evaluated  in  the  first-stage 
screen  using  just  a  limited  number  of  animals  (n  *  24);  a  similarly  small 
group  of  standard  decontaminant  animals  (n  =  24)  is  tested  concurrently  with 
the  experimental  decontaminant  animals.  It  was  recognized  by  MREF  personnel 
in  1985  that  a  considerable  data  base  would  be  amassed  over  time  pertaining  to 
the  lethality  rates  associated  with  the  standard  decontamination  procedure. 
This  historical  information  can  be  incorporated  into  the  efficacy  comparisons 
with  the  experimental  decontaminants  to  considerably  increase  the  statistical 
sensitivity  of  these  comparisons. 

Initial  statistical  methods  were  adopted  in  1985  to  make  use  of  the 
historical  information.  The  agent  dose  for  the  screening  test  was  set  on  the 
basis  of  an  extensive  LDj,  study  with  standard  decontaminant  animals;  the 
nominal  standard  lethality  rate  is  thus  50  percent.  The  lethality  rates 
observed  in  the  concurrent  standard  decontaminant  animals  were  compared  to 
this  nominal  50  percent  value.  If  the  concurrent  lethality  rate  differs  from 
the  nominal  by  more  than  three  standard  deviations,  then  the  concurrent  test 
is  considered  suspect  and  both  the  standard  and  experimental  decontaminant 
tests  are  repeated.  If  the  concurrent  lethality  rate  is  within  three  standard 
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deviations  of  50  percent,  the  experimental  decontaminant  lethality  rates  are 
compared  to  a  fixed  50  percent  standard.  If  the  concurrent  standard 
decontaminant  lethality  rate  is  between  two  and  three  standard  deviations 
above  the  nominal,  the  experimental  decontaminant  lethality  rate  is  also 
compared  to  the  concurrent  standard  decontaminant  lethality  rate.  The 
experimental  decontamination  procedure  is  declared  to  be  inferior  to  tie 
standard  procedure  if  the  observed  lethality  rate  among  the  experimental 
decontaminant  animals  is  statistically  significantly  greater  (P  *  0.05)  than 
0.5  (and  the  concurrent  standard  decontaminant  lethality  rate  in  those 
instances  when  it  is  compared  to  the  concurrent  standard  also). 

The  present  statistical  procedures  thus  compare  the  experimertal 
decontaminant  lethality  rate  either  to  the  concurrent  standard  decontaminant 
lethality  rate  or  to  a  fixed  50  percent  standard.  This  is  an  all  or  nothing 
procedure,  with  a  discontinuity  in  the  decision  point.  The  current  work  was 
undertaken  to  build  on  these  ideas  and  to  refine  their,  to  arrive  at  a  procedure 
which  allows  for  compromises  between  the  all  or  nothing  use  of  the  historical 
data.  Under  the  updated  procedure,  the  experimental  decontaminant  lethality 
rate  is  compared  to  a  weighted  average  of  the  concurrent  standard 
decontaminant  lethality  rate  and  the  historical  lethality  rate.  The  weight 
associated  with  the  concurrent  standard  lethality  rate  increases  as  the 
variability  in  the  concurrent  rate  decreases  (i.e.,  as  more  concurrent  animals 
are  tested),  as  the  observed  time  to  time  variability  in  the  historica 
standard  lethality  rates  increases,  and  as  the  agreement  between  the 
concurrent  rate  and  the  historical  rate  decreases. 

The  updated  procedure  is  based  on  the  assumption  that  the  historical 
standard  lethality  rates  for  individual  tests  are  completely  randomly 
distributed  about  an  overall  lethality  rate,  with  no  shifts  or  systematic 
drifts  in  the  rates.  Control  chart  methods  are  recommended  as  the 
surveillance  procedure  to  monitor  the  validity  of  this  assumption.  If  a  shift 
or  a  systematic  drift  over  time  is  detected  in  the  historical  standard 
decontaminant  lethality  rates,  the  data  from  the  far  past  should  be  excluded 
and  only  the  more  recent  past  data  should  be  used  to  form  a  historical  average 
standan  lethality  rate,  to  be  averaged  with  the  concurrent  standard  lethality 
results.  While  specific  procedures  for  the  elimination  of  data  from  the 
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historical  data  set  are  not  provided,  the  control  chart  procedures  described 
in  this  report  can  be  employed  to  determine  when  a  change  in  the  overall 
lethality  rate  has  occurred,  indicating  the  need  for  elimination  of  data  from 
the  historical  data  set. 

It  is  also  necessary  to  establish  appropriate  methods  for 
redetermining  an  agent  LDSI  dose  for  standard  decontaminant  protected  animals 
for  the  case  when  the  ongoing  surveillance  procedure  detects  a  significant 
shift  in  the  standard  decontaminant  lethality  rate,  signaling  a  corresponding 
significant  shift  in  the  LDSI  level.  The  LDSI  determination  procedure  in  the 
current  version  of  MREF  Protocol  21  (May  1985)  calls  for  using  a  minimum  of 
three  replicates  with  40  animals  per  replicate  to  determine  the  LDS(.  If 
insufficient  numbers  of  groups  are  obtained  with  observed  lethality  rates 
strictly  between  0  and  1,  then  additional  replicates  may  be  required.  Thus, 
200  or  more  animals  might  be  used  to  determine  an  LDSI  value  and  associated  , 
confidence  limits. 

Procedures  have  been  developed  in  conjunction  with  work  carried  out 
for  MREF  Tasks  85-18  and  87-34  to  estimate  the  LDSI  with  acceptable  precision 
based  on  many  fewer  animals.  To  accomplish  this,  it  is  necessary  that  the 
test  animals  be  distributed  among  appropriate  percentiles  of  the  dose-response 
distributions  for  that  agent  and  treatment  regimen.  These  dose-response 
distribution  percentiles  should  be  centered  around  the  true  LDSi  dose,  with 
sufficient  spread  that  the  dose-response  distribution  slope  may  be  determined. 
Since  the  dose-response  relationships  are  either  a  priori  unknown  or  just 
partially  known  based  on  historical  data,  the  allocation  of  animals  to  agent 
doses  is  made  in  a  stagewise,  adaptive  fashion  as  more  and  more  information 
about  the  current  dose-response  relation  becomes  available. 

The  methods  developed  here  should  be  considered  as  modifications, 
refinements,  and  improvements  to  the  methods  that  have  been  used  in 
conjunction  with  the  previous  MREF  Protocol  21  screening  program.  The  basic 
methodological  concepts  and  approaches  have  remained  unchanged. 

Statistical  problems  occurring  in  the  first-stage  screening  test  are 
addressed  in  this  report.  The  first  problem,  addressed  in  Section  2.0,  is  the 
development  of  a  test  procedure  for  comparing  each  experimental  decontaminant 
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with  the  standard  decontaminant.  Computer  programs  to  implement  the 
recommended  methods  have  been  developed  and  are  documented  in  Appendices  A 
and  B. 

The  second  problem,  addressed  in  Section  3.0,  is  the  development  of 
control  chart  procedures  for  monitoring  the  standard  deccntaminant  lethality 
rates  over  time.  Section  4.0  addresses  the  problem  of  comparing  the  standard 
decontaminant  lethality  rates  observed  in  replicate  subsets  of  the  concurrent 
test.  The  results  of  the  comparison  determine  whether  these  replicates  can  be 
pooled  to  arrive  at  an  overall  concurrent  lethality  rate,  and  if  not  which 
replicates  differ  from  the  others. 

Section  5.0  addresses  the  problem  of  ’•  determining  an  ID5#  dose 
level.  It  discusses  procedures  and  associated  r  jmputer  programs  to  carry  out 
stagewise,  adaptive  dose  allocation  designs  when  redetermining  the  LD5i.  It 
also  discusses  specialized  probit  analysis  model  fitting  procedures,  based  on 
nonlinear  regression  analysis,  for  updating  the  estimated  standard 
decontaminant  dose-response  distribution  following  each  stage  in  the  LDj, 
study. 

2.0  THE  TESTING  PROBLEM 

Each  time  a  first-stage  screening  test  is  performed  for  a  set  of 
experimental  decontairinants,  observed  lethality  rates  at  the  standard 
iecontaminant  105#  dose  are  obtained  for  the  standard  decontaminant  and  for 
•ach  of  the  experimental  deccntaminants,  based  on  a  limited  number  of 
laboratory  albino  rabbits  (nominally  24  animals  for  each  decontaminant).  A 
statistical  model  for  the  lethality  data  associated  with  the  standard  and 
experimental  decontaminants  is  oescribed  in  Section  2.1  and  the  testinq 
problem  is  stated  in  terms  of  the  parameters  of  this  model.  Recommended 
testing  procedures  are  developed  in  Section  2.2  and  are  characterized  in 
Section  2.3. 
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2.1  Statement  of  the  Testing  Problem 

Each  time  a  first-stage  screen  is  performed  for  a  set  of 
experimental  decontaminants,  a  limited  number  (nc)  of  animals  receive  the 
standard  decontaminant  treatment  and  a  nominal  LDU  dose  of  agent.  Also, 
limited  numbers  of  animals  receive  each  of  the  experimental  decontaminant 
treatments  and  the  same  nominal  10M  dose  of  agent  determined  for  the  standard 
decontaminant.  Each  experimental  decontaminant  is  compared  to  the  standard 
decontaminant  in  a  separate  statistical  test.  Thus,  without  loss  of 
generality,  assume  that  there  is  only  one  experimental  decontaminant  being 
tested;  denote  the  number  of  animals  receiving  this  experimental  decontaminant 
by  nt. 

The  number  of  lethalities  obtained  with  the  standard  decontaminant 
(xc)  is  assumed  to  have  a  binomial  distribution  with  ne  trials  and  success 
probabi 1 ity  pe,  where 


arcsin(lpe)  •  arcsin(ipt)  ♦  S.  (1) 

is  the  long-term  lethality  rate  for  the  standard  decontaminant  and  5  is  a 
random  (block)  effect  associated  with  this  particular  first-stage  screening 
test. 

The  number  of  lethalities  for  the  experimental  decontaminant  (xt)  is 
assumed  to  have  a  binomial  distribution  with  trials  and  success  probability 
pt,  where 

arcsin(Jpt)  •  arcsin(lut)  *  5.  (c) 

i»,  is  the  long-term  lethality  rate  for  the  experimental  decont.aminant  and  S'  is 
the  same  random  (block)  ?ffect  as  for  the  concurrent  lethality  rate,  pt. 

The  random  effect  term  <5  is  included  in  the  uiode I  to  account  for 
'hose  sources  of  experimental  variation  that  simultaneously  affect  the  true 
-♦hility  rates  for  all  the  standard  and  experimental  decontaminants  that  are 
'  *  “d  at  the  same  time.  It  is  assumed  that  5  is  distributed  as  a  mixture  of 

’  "ormal  distributions:  a  normal  <0,<7/^?)  distribution  with  probability  1-e 


and  a  normal (O.a^*)  distribution  with  probability  e.  Selecting  small  values 
of  (7£J  and  e  and  a  larger  value  of  05?  provides  a  model  that  results  in  small 
values  of  5  the  majority  of  the  time,  but  also  allows  for  an  occasionally 
large  random  effect.  This  reflects  the  situation  where,  on  most  occasions, 
the  true  lethality  rates  are  close  to  their  long  run  average  values  but  on 
infrequent  occasions  may  differ  considerably  from  those  long  run  values. 

The  arcsin-square  root  transformation  is  utilized  in  the  models  for 
pe  and  p%  in  anticipation  of  applying  the  same  variance-stabilizing 
transformation  to  the  observed  lethality  rates.  Let  rc  •  xe/ne  and  rt  ■  xt/nt 
denote  the  observed  lethality  rates  for  the  standard  and  experimental 
decontarainants,  respectively.  Also,  let  ae*  denote  0.25/ne  and  <rt*  denote 
0.25/nt.  For  the  purpose  of  deriving  a  test  statistic  for  comparing  the 
standard  and  experimental  decontaminants,  it  will  be  assumed  that  the 
conditional  distribution  of  arcsin(lr{)  given  5  is  approximately 
normal (arcsin(1pt),0.25/nc)  or  normal  (arcsin(!/»e)  +5, 0*) .  Similarly,  it  will 
be  assumed  that  the  conditional  distribution  of  arcsin(lrt)  given  6  is 
approximately  normal (arcs in(l pt) ,0.25/nt)  or  normal  (arcsin(lpt)  ♦$, 0*) . 

The  purpose  of  the  first-stage  screen  is  to  eliminate  those  and  only 
those  experimental  decontaminants  frcm  consideration  that  are  obviously 
inferior  to  the  standard  decontaminant.  The  problem  is  to  test  the  null 
hypothesis  H,:  /it  i  versus  the  alternative  Ht:  and  to  fail  the 

experimental  decontaminant  if  the  null  hypothesis  can  be  rejected  in  favor  of 
the  alternative  that  the  experimental  decontaminant  is  inferior. 

2.2  Recommended  Testing  Procedures 

Before  beginning  the  development  of  a  test  procedure,  we  must  first 
determine  the  criteria  that  the  test  must  satisfy.  Consider  the  criterion 
that  "the  test  must  have  a  (conditional)  significance  level  of  a  ■  0.05, 
conditioning  on  the  value  of  the  random  effect  5  at  the  time  the  screening 
test  is  performed".  If  the  significance  level  a  *  0.05  is  to  be  attained  for 
every  individual  realization  of  5,  then  5  is  being  treated  as  a  fixed  nuisance 
parameter  and  50  the  information  concerning  the  random  behavior  of  5  cannot  be 
iscd.  Lehmann  (1059)  indicates  that  the  standard  two-sample  test  is  best  in 
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this  situation.  However,  the  two-sample  binomial  test  does  not  make  use  of 
the  historical  information  available  concerning  the  fluctuations  of  5,  and 
therefore,  it  is  not  fully  efficient.  For  this  reason,  the  above  criterion 
was  not  pursued. 

Consider  the  alternative  criterion  that  "the  test  must  have  a 
(unconditional)  significance  level  of  a  *  0.05  where  the  random  nature  of  6  is 
factored  into  the  significance  level"  and  assume  that  the  null  hypothesis  will 
be  rejected  if  arcsin(!rt)  >  k,  where  k  is  a  critical  value  to  be  determined. 

This  formulation  adopts  a  random  effects  viewpoint.  Namely, 
hypothesis  tests  to  compare  concurrent  standard  and  test  decontamination 
procedures  will  be  carried  out  numerous  times  throughout  the  course  of  the 
screening  program.  The  random  (block)  effect  5  will  vary  across  these  tests 
according  to  a  probability  distribution,  which  can  be  estimated  based  on 
historical  data  and  which  is  discussed  below.  If  we  select  a  test  at  random 
from  the  "population"  of  such  tests,  we  wish  that  it  has  specified  Type  1 
error  level  a  (e.g.,  a  •  0.05).  If  5  varies  regularly  and  randomly  across 
tests,  then  this  formulation  permits  us  to  utilize  the  information  about  the 
standard  decontamination  procedure  lethality  rates  observed  in  the  previous 
tests  to  obtain  a  more  precise  estimate  of  the  current  standard 
decontamination  procedure  lethality  rate.  This,  in  turn,  results  in  increased 
sensitivity  of  the  current  test  of  hypothesis  to  compare  the  response  rates  of 
the  current  standard  and  experimental  decontamination  procedures,  relative  to 
what  would  be  obtained  if  the  historical  data  were  ignored. 

The  price  for  this  is  that  this  test  procedure  may  perform  more 
poorly  (e.g.,  in  terms  of  significance  levels  that  are  too  large  or  too  small) 
for  values  of  5  that  are  substantially  discrepant  from  past  values.  The 
variable  weighting  scheme  proposed  in  this  report  accounts  for  the  possibility 
of  occasional  values  of  5  that  are  somewhat  discrepant  from  historically 
observed  values. 

Since  the  distribution  of  rc  does  not  depend  on  (the  parameter 
about  which  v.e  are  making  an  inference),  the  principal  of  coi.ui  tional  i  ty 
implies  that  k  should  be  determined  from  the  conditional  null  disc,  -hut ion  of 
arcsin(lrt),  given  rc.  Denote  this  conditional  null  distribution  by 
F(arcsin(!  r.)  |  rr) .  Then  k  should  be  set  equal  to  FlH<I(arcsin(J  r.)  |  r,) ,  the 
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100 ( 1 -a) th  percentile  of  the  conditional  null  distribution.  Since  the 
conditional  significance  level  is  equal  to  a  for  each  value  of  re,  the  overall 
marginal  significance  level  is  exactly  a.  However,  the  conditional  (on  5) 
significance  level  is  a  function  of  5,  which  can  be  significantly  greater  than 
a  for  values  of  6  beyond  the  primary  support  of  the  assumed  probability 
distribution  of  5,  although  such  values  of  6  should  occur  only  very  rarely 
based  on  the  previous  history  of  the  screen. 

Consider  first  the  special  case  when  e  ■  0,  so  that  6  is 
normal  (0,a^2) .  Then,  under  the  null  hypothesis  that  *  /*c ,  F(arcsin(l  rt)  |re) 
is  approximately  normal  (/i|(j,2)  ,  where 


lit  3  arcsin(l/ie)  +  w  (arcsin(4re)  -  arcsin(!/ie)) 

3  (1  -  w)  arcsin(l/*e)  +  w  arcsin(lrc), 

(3) 

j,2  3  <rt2  +  w  a2, 

(4) 

w  3  erg,2  /  {ffg,2  *  6 r2) . 

(5) 

Equation  (3)  demonstrates  that  the  mean  itt  is  a  weighted  average  of  the 
transformed  long-term  lethality  rate  for  the  standard  decontaminant 
(arcsin(l/<c))  and  the  transformed  observed  lethality  rate  for  the  standard 
decontaminant  based  on  the  concurrent  control  animals.  Note  that  ttie  weight 
given  to  the  currency  observed  lethality  rate  increases  as  the  variability  of 
the  random  effect  5  (ofc2)  increases  and  as  the  variability  of  the  transformed 
observed  lethality  rate  for  the  standard  decontaminant  {a2)  decreases. 

However,  the  weight  does  not  depend  on  re. 

Thus  for  s  »  0,  the  critical  value  is  k  3  /i,  +  1.645  cr, .  The 
critical  region  an:sin(lr»)  >  k  employing  this  critical  value  satisfies  the 
criterion  stated  at  the  beginning  of  this  section.  However,  for  large  values 
of  \6\lofa,  the  conditional  (on  5)  significance  level  can  be  as  large  as  1. 

Such  values  of  the  ratio  are  highly  unlikely  under  the  assumption  that  e  3  0, 
yet  they  are  of  considerable  concern  in  the  development  of  the  hypothesis 
testing  procedure.  This  suggests  that  a  single  normality  assumption  for  6  may 
not  reflect  the  true  state  of  prior  feelings  about  the.  performance  of  the  test 
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system.  It  is  for  this  reason  that  we  have  chosen  to  assume  that  5  is 
distributed  as  a  mixture  of  two  normal  distributions  with  s.  >  0.  Such  a  model 
allows  for  the  occurrence,  on  infrequent  occasions,  of  more  extreme  values  of 
<5  than  would  be  predicted  by  a  simple  normal  distribution  model. 

Under  the  mixture  model,  the  random  variable  S  can  be  written  as 

6  »  (I-I)  Y,  +  I  Y1#  (6) 

where  I  has  a  binomial  distribution  with  1  trial  and  success  probability  e  and 
Y,  and  Yj  have  independent  norma) (0,<7£2)  and  normal (O,^2)  distributions. 

Then  the  conditional  distribution  of  I  given  rc  is  binomial  with  1  trial  and 
success  probability 

e*  -  cR  /  (1  -  e •+  eR),  (7) 

where 

*(arcsin(Jre),-arcsin(Juc)  ,<r&2+<re2) 

R  ■  - - -  (8) 

^(arcsin{lrc);arcsinO/ic),<r^J+iyc2)  . 

^(x;/i,a2)  is  the  normal  density  function  with  mean  ji  and  variance  a2  evaluated 
at  x.  Further,  F(arcsin(J  rt)  |  re)  is  a  mixture  of  two  normal  distributions:  a 
normal  (nt,at2)  distribution  with  probability  I  -  e*  and  a  normal  ,a2) 
distribution  with  probability  e',  where 


H j  3  arcsin0^e)  +  w-(  (arcs 

i n ( J rc )  -  arcsin(Vc)) 

(9) 

3  (1 

-  W:)  arcsin(lur)  +  w-  arcsin(lr, 

a2  3  j2  *  w 

i  ffc2. 

(10) 

and 

wi  3  ^5i2  /  W 

♦o- 

(11) 

The 

critical  value 

k  is  then  the  95th 

percenti le  of 

this  mixture 

distribution 

and 

the  recommended 

test  procedure  is 

to  reject  the 

null  hypothesis  Ha: 

Pt  ^ 

jit  in  favor  of 

the  alternative  H,: 

>  ,i*e .  i  f 

arcsin(Jr,i  > 

k. 
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The  critical  value  k  can  be  determined  in  an  iterative  fashion, 
starting  with  k,  *  p  +  1.645a,  where  p  *  (l-e')p9  +  e*/*i  and  3  (l-e*)cr,  + 
e*al.  The  ith  iterative  solution  (k5 )  can  be  expressed  as  a  function  of  the 
previous  solution  (k^)  as  follows: 

k;  -  kN1  +  (0.95  -  Ffkj.Jr,))  /  f(ki.1|re),  (12) 

where 

F(ki-x |  rc)  *  (1-e*)  *(ki-i .'/*«. O  +  e*  *(k,.l;^,a12) ,  (13) 

f(ki-ilre)  *  (i-**)  Kki.t.'A.'a*)  +  c*  (14) 

${x;p,o2)  is  the  cumulative  normal  distribution  function  with  mean  p  and 
variance  cr2  evaluated  at  x,  and  f(x;p,o2)  is  the  normal  density  function  with  , 
mean  p  and  variance  <j2  evaluated  at  x.  The  iterative  process  should  be 
continued  until  F(k,|re)  is  sufficiently  close  to  0.95. 

It  should  be  noted  that,  under  the  mixture  model,  the  relative 
weighting  of  the  concurrent  and  historical  standard  decontaminant  responses 
depends  on  the  value  of  rc.  The  farther  re  is  from  pc,  the  larger  is  R .  This 
implies  that  increasingly  more  weight  is  given  to  the  normal {px,a2)  , 
distribution,  which  in  turn  implies  that  increasingly  more  weight  is  given  to 
the  current  rc  in  the  determination  of  the  concurrent  standard  decontaminaht 
response  rate. 

Based  on  a  limited  examination  of  the  information  available  to 
support  a  selection  of  the  parameter*"  e,  a^,  and  agx,  it  is  recommended  that 
the  values  e  *  0.1,  3  0.1,  and  ogx  3  0.4  be  used  initially.  It  should  be 

noted  that  the  primary  reasons  for  the  selection  of  these  values  are  that  they 
are  consistent  with  the  historical  database  and  they  appear  to  provide  a  test' 
procedure  with  desirable  overall  properties  as  illustrated  in  Section  2.3. 
Further  work  must  be  performed  to  develop  a  procedure  for  the  selection  of 
test  procedure  parameters,  allowing  these  parameters  to  vary  with  the  agent 
and  test  system.  The  following  example  illustrates  the  use  of  the  procedure 
for  nt  =  nc  3  24  and  pc  3  0.5. 
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Suppose  we  have  nt  *  nc  *  24,  pe  3  0.5,  e  3  0.1,  3  0.1,  and 

crgi  3  0.4.  Then: 

arcsin0/ie)  3  t/4  3  0.7854 

*(arcsinOrc);0. 7854, 0.1704) 

R  3  - 

*(arcsin(Jre);0. 7854, 0.0204) 

t*  3  0.1  R/ (0.9  +  0.1  R) 

3  R/(9  +  R) 

w,  *  0.01/0.0204  =•  0.4898 

Wj  *  0.16/0.1704  =  0.9389 

p,  3  0.7854  +  0.4898  (arcsin(lrc)  -  0.7854) 

Pi  3  0.7854  +  0.9389  (arcsin(Jrc)  -  0.7854) 

a,2  3  0.0104  +  0.4898  (0.0104)  3  0.0155 
<7,  3  0.1246 

ffi2  3  0.0104  +  0.9389  (0.0104)  -  0.020197 
<7X  3  0.1421 

Values  of  R,  e*,  pQ,  pi,  k,  and  si n2(k)  are  listed  in  Table  2.2.1  for  various 
values  of  rc.  Either  of  the  las*1  two  columns  of  Table  2.2.1  can  be  used  to 
easily  carry  out  the  procedure  by  rejecting  the  null  hypothesis  Ha:  pt  <>  pc 
in  favor  of  the  alternative  Hr :  pi  >  p e  if  arcsin(lrt)  >  k  or  if  rt  >  sin2(k). 

2.3  Characterization  of  the  Testing  Procedures 

In  this  section,  we  characterize  and  compare  the  performance  of  the 
recommended  test  procedure  with  that  of  the  current  test  procedure  and  the 
standard  two-sample  binomial  test  procedure.  The  current  test  procedure 
involves  the  following  steps: 
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TABLE  2.2.1.  EXAMPLE  CALCULATIONS  FOR  THE  PROPOSED  TEST  PROCEDURE 
WITH  nt  =  n  =«  24,  u  =>  0.5,  e  3  0.1,  ok*  3  0.1, 

AND  <rft  3  0.4 


rc 

R 

e* 

to 

to 

k 

sin2(k) 

0.00 

206,026.54 

1.00 

0.40 

0.05 

0.28 

0.08 

0.10 

35.62 

0.80 

0.56 

0.35 

0.67 

0.33 

0.20 

3.22 

0.26 

0.63 

0.48 

0.82 

0.53 

0.30 

0.86 

0.09 

0.68 

0.59 

0.89 

0.60 

0.40 

0.43 

0.05 

0.74 

0.69 

0.94 

0.65. 

0.50 

0.35 

0.04 

0.79 

0.79 

0.99 

0.70 

0.60 

0.43 

0.05 

0.83 

0.88 

1.04 

0.75 

0.70 

0.86 

0.09 

0.89 

0.98 

1.11 

0.80 

0.80 

3.22 

0.26 

0.94 

1.09 

1.23 

0.89 

0.90 

35.62 

0.80 

1.01 

1.22 

1.44 

0.98 

1.00 

206,026.54 

1.00 

1.17 

1.52 

1.57 

1.00 
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(1)  Calculate  the  trigger  statistic 

Z  *  (re-/ie)/(/ie(l-/,e)/nc)1/2.  (15) 

(2)  If  Z  i  2,  calculate  the  test  statistic 

Ti  3  (rt-/‘c)/(/‘e(1"/‘c)/nt)1/2  (16) 

and  reject  the  null  hypothesis  H.:  pt  i  /*e  in  favor  of  the 
alternative  HjS  pt  >  pc  if  Tj  >  1.645. 

(3)  If  Z  >  2,  calculate  the  test  statistic 


rt(l-rt)  re(l*rc) 

— - - •+ - 


and  reject  the  null  hypothesis  H-:  /it  £  in  favor  of  the 
alternative  fy:  if  T2  >  1.645. 

The  standard  two-sample  binomial  test  procedure  is  based  on  the  test  statistic 
T2  defined  above  regardless  of  the  value  of  the  trigger  statistic  and  the  null 
hypothesis  H,:  5  /*c  is  rejected  in  favor  of  the  alternative  >  /ic  if 

T2  >  1.645. 

In  Tables  2.3.1  and  2.3.2,  the  three  test  procedures  are 
characterized  in  terms  of  the  probabi 1 ity  of  rejecting  the  experimental 
decontaminant  under  the  assumption  that  /*e  =  0.5.  In  Table  2.3.1,  the  tabled 
values  are  the  exact  conditional  probabilities  (x  1,000)  of  rejecting  the 
experimental  decontaminant  conditioned  on  the  value  of  the  lethality  rate  (pc) 
in  effect  at  the  time  of  the  test.  The  input  parameters  used  for  the 
recommended  procedure  are  o fa  =  0.1,  e  =  0.1,  =  0.4.  It  is  assumed  that 

n.  =  nc  =  24  and  /ie  =  0.5. 


TABLE  2.3.1  CONDITIONAL  PROBABILITIES  (x  1,000)  OF  REJECTING 
EXPERIMENTAL  DECONTAMINANT  FOR  SPECIFIC  VALUES 
OF  pt  AND  pe  ASSUMING  THAT  nt  ■  ne  ■  24  and 

Me  3  °*5 


‘Top  value  is  for  recommended  test  procedure  (cr*,  =  0.1,  e  =  0.1, 
h%-0.4) 

“Middle  value  is  for  current  test  procedure 

“Bottom  value  is  for  two-sample  binomial  test  procedure 


TABLE  2.3.2.  UNCONDITIONAL  PROBABILITIES  (x  1,000)  OF  REJECTING 
EXPERIMENTAL  DECONTAMINANT  FOR  SPECIFIC  VALUES  OF 
ftt  AND  VARIOUS  DISTRIBUTIONS  FOR  THE  RANDOM  EFFECT 
5  WITH  aSl  *  0.4  ASSUMING  THAT  nt  ■  ne  ■  24  and 
/»e  *  °-5 


e  -  0 

a5» 

t 

«  0.05 

<*5t 

e 

-  0.10 

e  -  0.20 

05* 

0 

0.05 

0.10 

0.20 

0.05 

0.10 

0.20 

0.05 

0.10 

0.20 

0.5 

042* 

051 

057 

074 

055 

059 

077 

066 

068 

082 

03  lb 

043 

064 

074 

.  047 

067 

C72 

049 

063 

068 

056c 

055 

053 

053 

055 

051 

054 

054 

054 

053 

0.6 

194 

195 

207 

222 

202 

212 

223 

208 

224 

233 

186 

197 

210 

192 

192 

207 

190 

189 

193 

182 

184 

179 

188 

182 

176 

181 

179 

168 

183 

181, 

0.7 

515 

515 

504 

483 

511 

505 

492 

508 

508 

495 

547 

510 

468 

385 

499 

451 

379 

471 

436 

359 

427 

424 

429 

425 

419 

427 

425 

413 

418 

418 

0.8 

852 

835 

823 

779 

P35 

815 

787 

829 

816 

784 

884 

821 

752 

612 

801 

732 

608 

765 

711 

593 

733 

725 

732 

720 

725 

722 

719 

712 

718. 

709 

0.9 

989 

984 

980 

966 

982 

978 

965 

979 

972 

965 

981 

955 

936 

837 

940 

918 

825 

896 

880 

809 

951 

944 

942 

924 

938 

937 

917 

921 

918 

911 

’Top  value  is  for  recommended  test  procedure  (ax,  3  0.1,  e  =  0.1 
<75,  =0.4) 

b Middle  value  is  for  current  test  procedure 

'Bottom  value  is  for  two-sample  binomial  test  procedure 


The  values  in  Table  2.3.1  for  *  0.5  are  the  conditional 
probabilities  that  an  experimental  decontaminant  with  the  same  lethality  rate 
as  the  standard  decontaminant  would  be  rejected  (false  rejection  rate)  in  a 
particular  screening  test  for  various  values  of  pe.  Note  that  the  two-sample 
binomial  test  procedure  has  a  relatively  constant  false  rejection  rate  for  all 
values  of  pe.  The  false  rejection  rate  of  the  current  test  procedure 
decreases  rapidly  as  oe  becomes  smaller  than  0.5  and  increases  rapidly  as  pc 
goes  from  0.5  to  approximately  0.7  and  then  decreases  as  pe  goes  from 
approximately  0.7  to  0.9.  The  false  rejection  rate  for  the  recommended  test 
procedure  follows  the  same  pattern  as  that  for  the  current  test  procedure  but 
with  less  rapid  increases  and  decreases.  The  remainder  of  Table  2.3.1 
illustrates  that  the  current  and  recommended  test  procedures  have  better 
conditional  power  relative  to  the  two-sample  binomial  test  procedure  for 
values  of  pc  greater  than  0.5  and  worse  power  for  values  of  pc  less  than  0.5. 

In  Table  2.3.2,  the  tabled  values  are  estimates  of  the  unconditional 
probabilities  (x  1,000)  of  rejecting  the  experimental  decontaminant  where  the 
random  effect  5  (and  therefore  pc)  is  allowed  to  vary  according  to  an  assumed 
probability  distribution.  Again,  the  input  parameters  used  for  the 
recommended  procedure  are  ag,  3  0.1,  e  3  0.1,  ogl  3  0.4.  It  is  assumed  that 
nt  =  nc  3  24,  fte  3  0.5  and  that  ergl  3  0.4  for  all  the  assumed  distributions 
for  <5. 

Each  value  in  Table  2.3.2  is  the  result  of  10,000  replications  of 
the  following  process.  Generate  a  5  value  from  the  mixture  of  two  normal 
distributions  defined  by  ag1  3  0.4  and  the  values  of  ag^  and  e  at  the  top  of 
the  column.  Generate  independent  binomial  test  results  for  the  standard  and 
experimental  decoritaminant  using  the  values  of  pc  and  pt  defined  by  the 
a»c  3  0.5,  <xt,  and  <5.  Record  the  result  of  each  of  the  three  test  procedures 
based  on  this  simulated  test  data. 

The  values  in  Table  2.3.2  for  3  0.5  are  the  probability  that  an 
experimental  decontaminant  with  the  same  lethality  rate  as  the  standard 
decontaminant  would  be  rejected  (false  rejection  rr.te)  for  various 
distributions  of  the  random  effect  5.  These  values  illustrate  that  the  false 
rejection  rate  is  reasonably  controlled  by  all  three  test  procedures  in  the 
neighborhood  of  the  assumed  distribution  for  5  (0^  =  0.1,  e  =  0.1,  ag x  3  0.4). 


The  remainder  of  Table  2.3.2  illustrates  that  the  recommended  test  procedures 
have  better  power  relative  to  the  current  test  procedure  and  the  two-sample 
binomial  test  procedure  ir.  the  neighborhood  of  the  assumed  distribution  for  5 
(oft  *  0.1,  e  ■  0.1,  afa  *  0.4).  j  This  increased  power  is  the  motivation  for 
the  development  of  the  recommended  test  procedure. 


3.0  THE  CONTROL  CHART  PROBLEM 


Over  time,  a  data  base  is  accumulated  for  the  standard 
decontaminant,  consisting  of  th$  observed  lethality  rates  from  the  individual 
screening  tests.  ^  statistical  model  for  the  lethality  data  associated  with 
the  standard  decontaminant  is  -described  in  Section  3.1  and  the  control  chart 
problem  is  stated  in  terms  of  the  parameters  of  this  model.  The  problem  of 
estimating  the  model  parameters  from  the  historical  database  is  discussed  in 
Section  3.2.  Recommended  contrdl  chart  procedures  are  developed  in 
Section  3.3  and  characterized  iri  Section  3.4. 


3.1  Statement  i of  the  Control  Chart  Problem 


Each  time  a  first-staae  screening  test  is  performed  for  a  set  of 
experimental  decontaminants,  a  limited  number  of  animals  receive  the  standard 
decontaminant  treatment  and  a  nominal  LDSi  dose  of  agent.  Let 

k  =  the  number  of  first  stage  screening  tests  in  the  historical 
database, 

n|  =  the  number  of  animals ! receiving  the  standard  decontaminant  during 
the  ith  screening  test,  and 


x,  =  the  number  of  lethalities  observed  with  the  standard  decontaminant 
during  the  ith  screening  test. 

X;  is  assumed  to  have  a  binomial  distribution  with  n,  trials  and  success 
probabi 1 i ty  p; ,  where 

arcsin(Jp;)  =  arcsin(l/*)  +<5,.  (18) 

ji  is  the  long-term  lethality  rate  for  the  standard  decontaminant,  and  5,  is  a 
random  effect  associated  with  the  ith  screening  test. 


c 
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The  random  effect  term  6 ;  is  included  in  the  model  to  account  for 
all  factors  that  randomly  affect  the  true  lethality  rate  for  the  standard 
decontaminant  from  test  to  test.  It  is  assumed  that  each  St  has  a 
Normal (O,^)  distribution  and  that  the  , * s  associated  with  separate  tests 
are  statistically  independent  of  one  another. 

The  arcsin-square  root  transformation  is  utilized  in  the  model  for 
p,  in  anticipation  of  applying  the  same  variance-stabilizing  transformation  to 
the  observed  lethality  rates.  Let  r,  •  Xj /n}  denote  the  observed  lethality 
rate  for  the  ith  screening  test.  Also,  let  a*  denote  0.25/n,.  For  the 
purpose  of  deriving  control  charting  procedures  for  monitoring  the  standard 
decontaminant  lethality  rate  over  time,  it  will  be  assumed  that  the 
conditional  distribution  of  arcsin(lrt),  given  5,  is  approximately 
normal  (arcsin(l  p;) , 0.25/n;)  or  normal  (arcsin(l^+5;  ,<y,) . 

The  purpose  of  the  control  chart  procedures  is  to  monitor  the 
standard  decontaminant  lethality  rate  over  time  to  detect  any  shifts  or  trends 
that  may  occur  as  a  result  of  random  and  inadvertent  variations  in  methods, 
materials,  or  agent  doses  employed.  The  problem  is  thus  to  plot  (a 
standardized  version  of)  rlf  rh  versus  time  along  with  upper  and  lower 

control  limits  that  characterize  the  expected  extreme  variatidns  according  to 
the  statistical  model  and  associated  parameter  estimates.  Values  beyond  the 
control  limits  are  evidence  that  either  the  statistical  model  or  the  parameter 
estimates  being  employed  may  no  longer  be  valid.  Aspects  of  the  tests  may 
then  need  to  be  adjusted,  for  example,  by  adjusting  the  agent  dose. 

3.2  Estimating  Model  Parameters 

The  first  step  in  forming  a  control  chart  for  the  standard 
decontaminant  legality  rate  is  to  determine  the  values  of  ft  and  that  will 
he  assumed  in  the  statistical  model.  The  following  procedure  may  be  used  to 
estimate  these  parameters  from  the  historical  database.  Let 

p  3  s i n2 (M) , 

E  w;  arcsin(lr) 

M  3 


(19) 


19 


and 


1 


o  i 


i 


(20) 


w;  Is  an  approximation  to  the  Inverse  of  the  variance  of  arcsln(lrl).  9/  Is 
defined  below  In  equation  (22).  Then  £  Is  an  unbiased  estimate  of  p  with 
approximate  standard  deviation 


Slip)  •  2  | 


sln(M)  cos (M) | 


(1/rw,)1'*. 


(21) 


Let 


where 


t  z,  9  * 
9#  -  - 

r  *1 


(22) 


(23) 


and 

9*  *  (arcsin(1r,)  -  H)*  -  (24) 


is  an  approximation  to  the  inverse  of  the  variance  of  9*.  Then  9$  is  an 
approximately  unbiased  estimate  of  of  with  approximate  standard  deviation 

SF {9j)  -  (2?) 


All  summations  above  are  over  i  •  l . k. 

Because  equations  (20),  (23),  and  (24)  involve  the  parameter 
estimates,  equations  (19)  and  (22)  must  be  solved  in  an  iterative  fashion. 
The  following  procedure  may  be  employed.  Begin  with  initial  values,  say 
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M  ■  ir/4  •  0.7854  (jj  ■  0.5)  and  •  0.01,  and  solve  equations  (19)  and  (22) 
using  these  values  in  equations  (20),  (23),  and  (24).  Repeatedly  solve 
equations  (19)  and  (22),  substituting  the  estimates  from  the  previous 
iteration  into  equations  (20),  (23),  and  (24)  until  the  estimates  do  not 
change  appreciably  from  one  iteration  to  the  next. 

3.3  Recommended  Control  Chart  Procedures 


As  stated  previously,  the  purpose  of  the  control  chart  procedures  is 
to  monitor  the  standard  decontaminant  lethality  rates  over  time  to  provide  a 
timely  signal  in  the  event  that  the  statistical  model  or  the  parameter 
estimates  being  employed  for  the  standard  decontaminant  are  no  longer  valid. 
Since  the  number  of  animals  used  may  vary  from  test  to  test,  it  is  most 
convenient  to  standardize  the  individual  transformed  lethality  rates  to 
achieve  approximate  uniform  variance  over  time.  Let 


Z. 


arcsInOr,)  -  arcsin(l^) 


0.25 


\  it 


(26) 


where  and  a?*  are  the  assumed  values  of  the  parameters  and  <j$.  Then 

Z,,...^  are  distributed  approximately  as  independent  standard  normal  random 

variables.  !t  is  recommended  that  lx . Zk  be  plotted  (as  the  vertical 

variable)  versus  time  (as  the  horizontal  variable)  along  with  horizontal  lines 
across  the  entire  plot  at  -3.00,  0,  and  3.00.  Three  types  of  control  chart 
procedures  are  considered: 


(A)  If  the  current  Z-value  falls  below  -3.00  or  above  3.00,  the  test 
involving  this  observation  should  be  repeated.  If  the  repeated 
test  also  exceeds  these  limits,  the  cause  of  this  exceedance 
should  be  investigated  and  corrected. 
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(B)  If  the  three  previous  Z-values  fall  either  all  below  -1.22  or 
all  above  1.22,  this  exceedance  is  strong  evidence  that  a  shift 
in  the  process  has  occurred.  The  cause  of  this  shift  should  be 
investigated  and  corrected. 

(C)  If  the  seven  previous  Z-values  fall  either  all  below  -0.28  or 
all  above  0.28,  this  exceedance  is  strong  evidence  that  a  shift 
in  the  process  has  occurred.  The  cause  of  this  shift  should  be 
investigated  and  corrected. 

These  three  procedures  provide  short-term,  intermediate,  and 
long-term  tests  respectively  for  a  shift  in  process  behavior.  As  with  the 
testing  procedure,  it  is  recommended  that  the  value  og,  be  set  to 
0.1  initially. 

If  any  of  the  three  control  chart  procedures  exceed  thei”  critical 
values,  the  historical  database  should  be  scrutinized.  If  some  of  the  older 
data  are  no  longer  pertinent  to  current  tests,  they  might  be  eliminated  from 
calculations  of  model  parameters.  Consideration  might  be  given  to  carrying 
out  a  new  LDSI  study  and  adjusting  the  agent  dose. 

The  above  procedures  are  designed  to  be  easy  to  carry  out  with  a 
calculator  and  a  simple  plot  of  the  data.  Similar  tests  based  on  the  median 
can  be  employed  as  follows.  Let 


Mt  ■  the  Z-value  (Equation  26)  for  the  most  previous  test 

Mj  ■  the  median  of  the  Z-values  (Equation  26)  for  the  three 
previous  tests,  and 

W,  *  the  median  of  the  Z-values  (Equation  26)  for  the  seven 
previous  tests. 


The  following  procedures  are  analogous  to  procedures  A,  B,  and  C  defined 
above. 
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(A2)  If  Mj  falls  below  -3.00  or  above  3.00,  the  test  involving  this 
observation  should  be  repeated.  If  the  repeated  test  also 
exceeds  these  limits,  the  cause  of  this  exceedance  should  be 
investigated  and  corrected. 

(B2)  If  M3  falls  below  -2.17  or  above  2.17,  this  exceedance  is 

strong  evidence  that  a  shift  in  the  process  has  occurred.  The 
cause  of  this  shift  should  be  investigated  and  corrected^ 

(C2)  If  falls  below  -1.42  or  above  1.42,  this  exceedance  is 

strong  evidence  that  a  shift  in  the  process  has  occurred.  The 
cause  of  this  shift  should  be  investigated  and  corrected. 


Similar  tests  based  on  the  mean  can  be  employed  as  follows.  Let 

Tj  *  the  Z-value  (Equation  26)  for  the  most  previous  test 

Tj  *  the  average  of  the  Z-values  (Equation  26)  for  the  three 
previous  tests,  and 

Tr  *  the  average  of  the  Z-values  (Equation  26)  for  the  seven 
previous  tests. 


The  following  procedures  are  analogous  to  procedures  A,  B,  and  C  and  A2,  B2, 
and  C2  defined  above. 


(A3)  If  Tt  falls  below  -3.00  or  above  3.00,  the  test  involving  this 
observation  should  be  repeated.  If  the  repeated  test  also 
exceeds  these  limits,  the  cause  of  this  exceedance  should  be 
investigated  and  corrected. 

(B3)  If  T3  falls  below  -1.73  or  above  1.73,  this  exceedance  is 

strong  evidence  that  a  shift  in  the  process  has  occurred.  The 
cause  of  this  shift  should  be  investigated  and  corrected. 

(C3)  If  T7  falls  below  -1.13  or  above  1.13,  this  exceedance  is 

strong  evidence  that  a  shift  in  the  process  has  occurred.  The 
cause  of  this  shift  should  be  investigated  and  corrected. 
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3.4  Characterization  of  the  Control  Chart  Procedures 


In  Table  3.4.1,  the  three  sets  of  control  chart  procedures  [ (A , 8 , C) , 
(A2,B2,C2),  and  (A3 , B3 , C3) 3  are  characterized  in  terms  of  the  probability  of 
signaling  a  process  shift.  The  model  parameters  employed  are  fie9  *  0.5  and 
fffa  •  0.1  and  it  is  assumed  that  n;  *  24. 

The  values  in  Table  3.4.1  are  the  probability  that  a  process  shift 
will  be  detected  by  t  ?  various  procedures.  The  results  are  based  on  the 

assumption  that  the  Z  values  are  independent  standard  normal  random  variables. 

The  results  for  the  B l  and  C2  procedures  are  based  on  a  normal  approximation 

to  the  distribution  or‘  median  of  independent  standard  normal  random  variables. 

While  the  procedures  based  on  counts  and  medians  are  simpler  to  carry  out,  it 
is  recommended  that  the  procedures  based  on  means  (A3,  B3,  and  C3)  be  employed 
due  to  the  significant  power  advantage  for  detecting  moderate  shifts  in  the 
standard  decontaminant  lethality  rate. 

4.0  COMPARISONS  OF  THE  STANDARD  DECONTAMINANT  LETHALITY  RATES 
AMONG  REPII CATE$ . WITHIN  "TESTS' 

Each  screening  test  involves  the  simultaneous  testing  of  nc  animals 
with  the  standard  decontaminant  and  nt  animals  with  each  of  the  experimental 
decontaminants.  Usually  nc  -  »  24.  For  logistical  reasons,  particularly 

if  a  number  of  test  decontaminants  are  to  be  evaluated  at  the  same  time,  the 
test  is  divided  into  K  replicate  portions  and  each  portion  is  carried  out  on 
separate  days.  Usually  K  *  3,  and  n  *  8  animals  are  tested  per  group  per  day. 

Comparisons  between  the  standard  decontaminant  results  and  the  test 
decontaminant  results  usually  incorporate  the  assumption  that  the  individual 
replicate  results  within  tests  can  be  pooled  to  arrive  at  overall  lethality 
rates.  Preliminary  comparisons  among  the  standard  decontaminant  lethality 
rates  observed  in  each  replicate  are  carried  out  to  examine  the  reasonableness 
of  this  assumption.  If  there  is  no  evidence  of  heterogeneity  among  the 
standard  decontami nant  replicates,  then  it  is  presumed  that  the  replicates 
■vere  carried  out  under  homogeneous  conditions  and  the  test  decontaminant 
results,  as  well  as  the  standard  decontaminant  results,  are  pooled  across 


s 
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TABLE  3.4.1.  PROBABILITIES  (x  1,000)  OF  SIGNALING  A  PROCESS 
SHIFT  FOP  SPECIFIC  VALUES  OF  ne  ASSUMING  THAT 
rij  *  24,  net  3  0.5,  AND  Ofa  3  0.1 


Control  Chart 

Procedure 

Me 

A 

B 

C 

0.5 

003* 

003 

003 

003b 

003 

003 

003c 

003 

003 

0.6 

Oil 

028 

057 

on 

021 

066 

Oil 

038 

128 

0.7 

059 

202 

399 

059 

157 

517 

059 

306 

791 

0.8 

227 

612 

842 

227 

545 

960 

227 

816 

998 

0.9 

597 

937 

989 

597 

931 

1,000 

597 

996 

1,000 

’Top  value  is  for  procedures  A,  B,  and  C  (Counts) 
bMiddle  value  is  for  procedures  A2,  B2,  and  C2  (Medians) 
'Bottom  value  is  for  procedures  A3,  B3,  and  C3  (Means) 
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replicates.  The  pooled  standard  decontaminant  results  are  then  compared  with 
the  historical  standard  decontaminant  results  and  with  the  control  chart 
limits,  in  the  manner  discussed  in  Sections  2.0  and  3.0.  The  resulting 
estimate  of  current  standard  decontaminant  lethality  rate  is  compared  with 
each  of  t;he  test  decontaminant  rates,  in  the  manner  discussed  in  Section  2.2. 

If  there  is  statistically  significant  heterogeneity  among  the 
standard  decontaminant  replicates,  further  tests  are  carried  out  to  determine 
which  replicates  differ  from  the  others  and/or  from  the  long  run  historical 
results.  These  outlying  replicates  are  then  considered  for  deletion  (both  the 
standard  decontaminant  and  the  test  decontaminant  results)  and  additional 
replicates  are  carried  out  to  replace  them. 

4.1  Comparison  of  Individual  Replicate  Standard  Decontaminant 
~  Resuits  with  the  Overall  Average  ~T"' 


Suppose  that  the  current  test  is  divided  into  K  replicates  (days), 
that  the  ith  replicate  includes  n;  animals  in  the  standard  decontaminant 
group,  and  that  x,  responses  (deaths)  are  observed  among  these  animals.  Let 
r j  =  Xj/n;  denote  the  observed  response  rate  in  the  ith  replicate,  xc  *  1^,  jX;, 
nc  =  n  =»  nc/K,  and  rc  =  xe/ne.  Let  Pj  *  E( r; )  denote  the  population 

average  response  rate  in  the  ith  replicate.  The  analysis  of  means  (Ott,  1975) 
is  used  to  compare  each  replicate  response  rate,  rjt  to  the  average  rate,  rc. 
The  analysis  of  means  test  is  designed  to  be  sensitive  to  the  presence  of  an 
extreme  replicate  that  differs  from  the  others,  much  like  a  control  chart 
inference.  The  hypothesis 

H0:  Pi  3  P2  3  •  •  •  "  Pk 
is  tested  by  the  analysis  of  means  procedure.  Let 

Z  =  max  [| r;  -  rc|/(rc(l  -  rc)/n)*] 
i  3  1 . K 
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The  hypothesis  H0  is  rejected  at  significance  level  sc  if  Z  >  fy,  where  Hj  is 
tabulated  by  Ott  (1975),  Schilling  (1973),  and  others.  If  K  =*  3  and  a  *  0.05, 
then  Ha ls  *  1.93. 

If  H0  is  rejected,  each  rj  is  compared  to  the  others  and  to  the 
historical  control  rate  (nominally  0.50) . 

4.2.  Comparisons  of  Replicate  Responses  with  the  Historical 
Control  Rate  and  Among  Each  Other' 

If  the  analysis  of  means  test  rejects  the  above  H0,  then  each  r,  is 
compared  to  the  long  run  historical  standard  decontaminant  response  to 
determine  which  differ.  Assume  for  purposes  of  this  discussion  that  the 
historical  response  rate  is  close  to  the  nominal,  0.50.  The  hypotheses 

Hb:  pj  =*  0.50  and  i  «  1,  2,...,K 

are  tested  using  individual  one-sample,  two-sided  tests.  The  significance 
levels  are  adjusted  by  Bonferroni's  method  so  that  the  overall  Type  1  error 
across  the  K  tests  does  not  exceed  cc.  Let 

Zj  ■  \r,  -  0.5|/[(0.5)(0.5)/n,]* 

The  hypothesis 

H,:  Pi  -  0.5 

is  rejected  if  Zt  >  Zg.  Values  of  for  a  =  0.10,  0.05,  and  0.01  and 

K  =  1(1)  10  and  above  are  tabulated  by  Ott  (1975),  Miller  (1966),  and  others. 

If  K  =  3,  then  Z,  ,5  =  2.39  and  Z,  lg  =  2.11. 

It  should  be  noted  that  if  K  =  3  and  n,  =  8,  only  r-  -  0  or  r-  =  1 
would  cause  H0  to  be  rejected  at  a  =  0.05;  r,  l  0.875  or  r;  i  0.125  would 

cause  H0  to  be  rejected  at  a  3  0.10.  When  ^  =  8,  these  tests  are  rather 

insensitive;  they  will  detect  only  very  large  departures  from  consistency 
across  replicates  or  from  the  historical  average  rate. 


27 


Pairwise  comparisons  among  replicates  are  carried  out  using  Fisher's 
exact  test  (two-sided).  With  just  n;  -  8  animals  per  group,  these  tests  are 
also  insensitive.  Critical  values  are  tabulated  by  Pearson  and  Hartley  (1958) 
When  K  ■  3,  the  single-tailed  0.01  level  critical  values  correspond  to  an 
a  *  0.06  (0.01  x  2  x  3)  two-tailed  simultaneous  significance  level  for  all 
(three)  pairwise  comparisons  among  replicates.  The  Pearson-Hartley  table 
demonstrates  that  8  of  8  responses  can  be  distinguished  from  2  of  8,  7  of  8 
from  1  of  8,  and  6  of  8  from  0  of  8  at  this  significance  level.  The  single¬ 
tailed  0.025  level  critical  values  correspond  to  an  a  *  0.12  (0.020  x  2  x  3) 
two-tailed  simultaneous  significance  level.  At  this  significance  level, 

8  of  8  can  be  distinguished  from  3  of  8,  7  of  8  from  2  of  8,  6  of  8  from 
1  of  8,  and  5  of  8  from  0  of  8.  Thus,  when  n;  *  8,  these  pairwise  comparisons 
will  detect  omy  substantial  departures  from  consistency  across  replicates. 

The  discussion  in  this  section  demonstrates  that  only  sizeable 
inconsistencies  among  replicate  response  levels  will  be  flagged  by  these 
procjdures.  In  all  other  instances,  the  responses  will  be  pooled  across 
replicates  and  the  principal  comparisons  will  proceed. 

5.0  REDETERMINATION  OF  LD„  DOSES 

An  important  aspect  of  the  screening  program  is  to  initially 
establish  and  then  periodically  update  L058  doses  for  standard  treatment.  If 
the  control  chart  inferences  discussed  in  Section  3.0  detect  drift  in  the 
standard  decontaminant  response  rates  or  repeated  exceedences  of  the  control 
chart  limits,  then  the  LDsg  dose  needs  to  be  redetermined.  A  new  LD5a  study 
must  be  carried  out  to  determine  the  new  agent  dose.  This  section  discusses 
experimental  design  and  data  analysis  methods  and  associated  computer  programs 
that  have  been  developed  to  determine  LDsa  doses  in  an  efficient  manner;  fewer 
animals  are  needed  to  attain  the  desired  levels  of  estimation  precision, 
relative  to  a  classical  LD5B  design. 

Section  5.1  discusses  a  stagewise  dose  allocation  experimental 
design  strategy  that  has  been  developed  to  accomplish  this  aim.  Such 
stagewise  designs  lead  to  nontraditional  dose  allocations  that  utilize 
relatively  large  numbers  of  doses  with  relatively  small  numbers  of  animals  per 
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dose.  It  is  possible  that  each  animal  will  be  tested  at  a  different  dose. 
Standard  probit  analysis  computer  programs,  therefore,  cannot  be  used  to  fit 
dose-response  models  to  the  lethality  data. 

Specialized  procedures,  based  on  nonlinear  regression  analysis,  have 
been  developed  to  fit  dose-response  models  to  these  data,  "hese  procedures 
have  been  developed  in  a  series  of  computer  programs  based  on  the  general 
purpose  nonlinear  regression  procedure,  PROC  NUN,  in  the  SAS  statistical 
computing  system  (SAS,  1985).  Section  5.2  discusses  the  procedures  and 
programs. 

5.1  Staqewise,  Adaotive  Dose  Allocation  Procedures 

The  LDSi  is  estimated  based  on  a  small  to  moderate  number  of 
animals.  The  precision  of  estimation  of  the  LDS>  and  slope  depends  on  the 
numbers  of  animals  tested  as,  well  as  on  the  allocation  of  animals  to 
appropriate  portions  of  the  (unknown)  dose-response  distribution.  To  obtain 
relatively  precise  estimates  of  the  LD5I  and  slope  with  the  numbers  of  animals 
available,  the  test  doses  should  be  centered  around  the  LDS,  with  enough 
spread  to  permit  good  estimation  of  the  slope.  The  test  doses  should  not, 
however,  extend  too  far  beyond  the  central  portion  of  the  dose-response  region 
(e.g.,  they  should  lie  between  the  10th  and  90th  percentiles).  The  desired 
dose  allocation  heavily  depends  on  the  underlying  dose-resppnse  distribution. 
It  is  assumed  that  the  dose-response  relation  for  the  standard  decontaminant 
animals  can  be  described  by  a  two-parameter  probit  model  without  background, 
at  least  in  the  central  portion  of  the  dose-response  region. 

The  relative  sensitivities  of  alternative  dose  allocations  can  be 
evaluated  before  any  data  have  been  collected.  This  permits  "target  designs" 
to  be  selected  before  the  start  of  the  experiment  and  to  be  updated  as  the 
experiment  proceeds. 

Since  the  underlying  dose-response  distributions  are  not  known  prior 
to  the  start  of  the  test,  the  LDSI  test  is  carried  out  in  a  stagewise  fashion. 
The  dose  allocation  for  the  first  stage  is  based  on  historical  results. 
Previous  LD5g  study  results,  augmented  by  observed  response  rates  in  more 
recent  standard  decontaminant  tests.,  can  be  used  to  obtain  initial  estimates 
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of  the  LDSf  and  slope  and  the  associated  first-stage  dose  allocation. 

Following  each  stage,  the  information  concerning  the  underlying  dose-response 
distribution  is  updated  based  on  fitting  dose-response  models  to  the  results 
obtained  in  the  current  and  previous  stages.  Doses  are  selected  for  the  next 
stage  to  best  approximate  the  target  design,  over  and  above  the  previous 
allocations,  based  on  the  updated  dose-response  distribution.  This  process  is 
iterated  until  the  completion  of  all  the  stages  or  until  the  LDS#  and/or  slope 
are  estimated  with  the  required  level  of  precision. 

This  approach  is  in  the  spirit  of,  but  is  more  flexible  and  adaptive 
than,  the  formal  up-dcwn  method  (Dixon  and  Mood,  1948).  It  attempts  to 
incorporate  relatively  large  numbers  of  test  doses  within  each  stage  and  uses 
information  from  all  previous  stages  to  make  decisions  about  the  doses  to  be 
selected  in  suosequent  stages. 

The  application  of  the  stagewise  dose  allocation  approach  to 
determining  the  standard  decontaminant  LDSI  is  illustrated  by  an  example 
pertaining  to  percutaneous  application  of  GD  in  albino  rabbits  and  treatment 
with  both  components  of  the  M258A1  standard  decontamination  kit-  An  updated 
program  to  screen  new  candidate  decontaminants  is  to  be  implemented,  utilizing 
whatever  information  can  be  obtained  from  previously  completed  screening 
programs. 

The  a  priori  assumptions  for  the  dose-response  relation  applicable 
to  the  forthcoming  screening  program  are  based  on  a  probit  model  fit  to  the 
results  from  a  previous  LDSi  study  with  this  same  agent,  animal  model,  and 
decontamination  regimen  that  was  carried  out  in  May-June  1985  in  and  MREF 
Final  Report  entitled  "Task  85-10:  Validation  of  a  Protocol  to  Compare  the 
Effectiveness  of  Experimental  Decontaminants  With  Both  Components  of  the 
M258A1  Kit  Against  Percutaneous  Application  of  Undiluted  Organophosphate 
Chemical  Surety  Materiels  t.o  the  Laboratory  Albino  Rabbit,"  (December  1987, 
Table  3.1.4).  The  LDSJ  was  estimated  to  be  13.0  ^g/kg  and  the  slope  3.732, 
based  on  n  =  360  animals.  The  parameters  of  this  distribution  are  displayed 
in  Table  5.1.1. 

Ten  alternative  "target  designs"  were  considered.  These  are' 
numbered  GDI  to  GD10  and  are  shown  in  Table  5.1.2.  Each  target  design 
consists  of  n  =  100  animals,  allocated  equally  or  unequally  to  various 


TABLE  5.1.1.  PARAMETERS  SPECIFYING  THE  A  PRIORI  PROBIT  DOSE-RESPONSE 
DISTRIBUTION.  THESE  PARAMETERS  DETERMINE  THE  CENTRAL 
DISTRIBUTION  AND  PERTURBATIONS  INCORPORATED  IN  SUBSEQUENT 
SENSITIVITY  ANALYSES 


FITID  BO  B1  VO  COl  VI 

GD/OECON  0,843  3.732  0.336  -0.299  0.274 


FITID  =•  Fit  identification 

BO , B1  =  Slope  and  intercept  of  the  a  priori  dose-response  distribution 
VO, VI, COl  *  Variance  of  the  intercept,  variance  of  slope,  and  covariance 
between  the  intercept  and  slope  of  the  a  priori  dose-response 
distribution  and  quantify  the  uncertainty  in  these  parameters. 
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combinations  of  the  10,  20,... ,90  percentiles  of  the  assumed  prior 
dose-response  distribution.  For  example,  target  design  GDI  allocates 
25  animals  to  each  of  the  20,  *0,  60,  and  80  percentiles.  The  absolute 
sensitivities  calculated  for  these  designs  pertain  to  100  animals.  However, 
these  sensitivities  are  each  scaled  up  or  down  by  the  factor  (100/n)^  if  n 
animals  are  used  instead.  The  rel atl ve  sensitivities  for  these  designs  are 
thus  invariant  to  sample  size. 

Although  each  target  design  is  presented  for  initial  planning 
purposes  as  a  single-stage  design,  the  LDM  study  is  in  fact  carried  out  in 
stages,  with  doses  adjusted  from  stage  to  stage.  The  target  allocations  are 
updated,  over  and  above  the  doses  previously  tested,  in  light  of  the  most 
current  estimate  of  the  dose-response  relation.  The  stagewise,  adaptive  dose 
allocation  helps  assure  conformance  to  the  target  design  even  if  the  estimate 
of  the  underlying  dose-response  distribution  shifts  from  stage  to  stage  as 
additional  results  are  obtained.  Furthermore,  if  the  attained  sensitivity  to 
estimate  dose-response  distribution  parameters  or  to  compare  dose-response 
distributions  exceeds  that  predicted  at  the  outset  of  the  experiment,  the 
stagewise  design  strategy  can  lead  to  early  stopping. 

The  predicted  sensitivities  for  each  design  are  calculated  from  the 
information  obtained  in  the  previous  stages,  combined  with  the  expected 
information  to  be  obtained  in  the  current  and  future  stages.  The  information 
associated  with  each  design  is  evaluated  for  the  distribution  specified  in 
Table  5.1.1  (the  "central"  distribution),  as  well  as  tor  distributions  that 
are  perturbations  about  the  central  distribution. 

Table  5.1.3  displays  the  "results"  from  the  "previous"  stages.  For 
eacn  stage  and  dose,  the  logarithm  of  dose  (X),  the  number  of  animals  on  test 
(Nfl) ,  and  the  number  of  responses  (Y)  are  given.  In  this  example,  the 
designs  are  being  evaluated  prior  to  the  first  stage.  There  is  no  previous 
data  and  so  NN  *  0.  For  evaluations  following  later  stages,  the  observed 
results  at  all  the  previous  stages  and  doses  would  be  used. 

Detailed  sensitivity  analyses  are  carried  out  for  each  target  design 
to  assess  its  performance  under  a  variety  of  distributions  that  niaht  be 
lively  to  occur  (i.e.,  perturbations  about  the  central  distribution). 
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TABLE  5.1.3.  EXPERIMENTAL  "RESULTS*  FROM  "PREVIOUS"  STAGES1 


OBS 

GROUP 

STAGE 

OOSE 

X 

NN 

Y 

1 

DECON 

1 

20 

1.30103 

0 

0 

•Prior  to  the  first  stage,  the  numbers  of  animals  (NN)  and  the  numbers  of 
responses  (Y)  are  each  0. 

OBS  «  Record  number 

GROUP  -  Identification  variable 

STAGE  •  Stage  at  which  dose  was  administered 

DOSE  *  Oose  administered 

X  »  Common  logarithm  of  dose 
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Table  5.1.4  displays  the  detailed  sensitivity  results  for  design  GDI.  The 
middle  line  in  the  table  (underlined)  corresponds  to  the  central  distribution. 
The  rt  lining  48  lines  correspond  to  perturbations  about  this  central 
distribution.  Standard  errors  of  the  estimates  of  the  specified  logarithmic 
percentiles  (50,  80,  90  in  this  example)  and  of  the  slope  are  calculated  for 
each  of  these  distributions.  The  results  of  the  sensitivity  analyses  for  each 
target  design  are  summarized  in  Tables  5.1.5  and  5.1.6.  Table  5.1.5  displays 
weighted  averages  of  the  standard  errors  over  all  of  the  49  distributions  in 
the  sensitivity  analysis;  the  distributions  closer  to  the  central  distribution 
receive  the  greater  weight.  Table  5.1.6  displays  the  minima  and  the  maxima  of 
these  standard  errors  over  these  same  distributions.  The  maxima  can  be 
regarded  as  "worst  cases",  over  the  range  of  distributions  considered 
plausible  based  on  the  current  information. 

Table  5.1.5  shows  that  the  weighted  averages  of  the  standard  errors 
of  the  log,,  (LD5,)  are  similar  across  all  of  the  target  allocations 
considered.  This  is  not  surprising,  since  they  were  all  selected  to  be 
symmetric  about  the  a  priori  50th  percentile.,  Design  GD6  has  the  smallest  and 
design  GD9  has  the  largest.  By  contrast  the  standard  errors  of  the  log,, 

(ID„)  and  the  slope  vary  to  a  greater  extent  across  the  target  allocations 
considered.  Oesign  GD9  has  the  smallest  and  design  GD6  is  the  largest. 

Similar  considerations  hold  for  the  maxima  of  the  standard  errors. 
Those  for  the  log,,  (LDS,)  are  similar  across  all  the  target  allocations 
considered.  Those  for  the  log,,  (ID„)  and  the  slope  vary  more  across  the 
target  allocations  considered;  design  GD9  has  the  smallest  and  design  GC6  has 
the  largest. 

Oesign  GD6  allocates  animals  evenly  to  the  30,  50,  and  70 
percentiles;  design  G09  allocates  animals  evenly  to  10,  50  and  90  percentiles. 
Designs  GD3  and  GD10  are  compromises  between  the  two  extremes.  They  allocate 
animals  equally  among  the  10,  30,  50,  70  and  90  percentiles  and  among  the  10, 
20,  50,  80,  and  90  percentiles,  respectively. 

To  utilize  this  information  for  a  stagewise  dose  allocation,  the 
numbers  of  stages  and  the  numbers  of  animals  per  stage  would  be  decided  upon 
and  design  G03,  for  examole,  appropriately  scaled  down,  might  be  run  for  the 
first  stage.  Following  the  first-stage,  the  dose-resDonse  distribution 


lAiilt  5.1.4.  Ct  FAILED  SENSITIVITY  ANALYSIS  RESULTS  TOR  DESIGN  GDI.  UNDERLINED 
DISTRIBUTION  IS  IDE  CENTRAL  DISTRIBUTION 
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TABLE  5.1.5.  WEIGHTED  AVERAGES  OF  THE  STANDARD  ERRORS  OF  THE  LOGARITHMIC 
50,  80,  AND  90  PERCENTILES  AND  THE  SLOPE  OVER  ALL  THE 
DISTRIBUTIONS  IN  THE  SENSITIVITY  ANALYSIS 


Weighted  averages  of  standard  errors  of  points  varied  around 
A  and  B  over  the  designs  in  the  sensitivity  analysis 


Design  I.D.  3  GDI 


APERCENT 

IDPCT 

NPCT 

SQRTSUM 

50 

,  SEFVX50 

49 

0.036863 

80 

SEFVX80 

49 

0.065796 

90 

SEFVX90 

49 

0.090660 

B1 

SEFVX31 

49 

0.850717 

Design 

I.D.  3  GD10 

50 

SEFVX50 

49 

0.038451 

80 

SEFVX80 

49 

0.058843 

90 

SEFVX90 

49 

0.077879 

B1 

SEFVXB1 

49 

0.704621 

Design 

1.0.  3  GD2 

50 

SEFVX50 

49 

0.037861 

80 

SEFVX80 

49 

0.061135 

90 

SEFVX90 

49 

0.082221 

B1 

SEFVXB1 

49 

0.759807 

Oesign 

I.D.  3  GD3 

50 

SEFVX50 

49 

0.038780 

80 

SEFVX80 

49 

0.058818 

90 

SEFVX90 

49 

0.077645 

B1 

SEFVXB1 

49 

0.705694 

* 

Design 

I.D.  =  GD4 

50 

SEFVX50 

49 

0.037343 

80 

SEFVX80 

49 

0.062219 

90 

SEFVX90 

49 

0.084367 

81 

SEFVXB1 

49 

0.778833 
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TABLE  5.1.5. 
(Continued) 


Design  I .0.  3  GD5 

■ 

APERCENT 

IDPCT 

NPCT 

SQRTSUM 

50 

SEFVX50 

49 

0.037063 

80 

SEFVX80 

49 

0.064222 

90 

SEFVX90 

49 

0.087923 

B1 

SEFVXB1 

49 

0.820809 

Design  I.D.  3  GD6 

50 

SEFVX50 

49 

0.03609 

80 

SEFVX80 

49 

0.08256 

,  90 

SEFVX90 

49 

0.11852 

B1 

SEFVXB1 

49 

1.14917 

Design  I.D.  3  GD7 

50 

SEFVX50 

49 

0.039233 

80 

SEFVX80 

49 

0.058390 

90 

SEFVX90 

49 

0.076798 

81  . 

SEFVX81 

49 

0.693325 

Design  I.D.  3  GD8 

50 

SEFVX50 

49 

0.038148 

80 

SEFVX80 

49 

0.059842 

90 

SEFVX90 

49 

0.079798 

B1 

SEFVX81 

49 

0.728632 

Design  I.D.  3  GD9 

50 

SEFVX50 

49 

0.040943 

80 

SEFVX80 

49 

0.056355 

90 

SEFVX90 

49 

0.071733 

B1 

SEFVXB1 

49 

.  0.623474 

APERCENT  3  Quantity  being  estimated  (percentile  or  slope) 

IDPCT  =  Identifier  of  quantity  being  estimated 

NPCT  =  Number  of  distributions  entering  into  the  weighted  average 

SQRTSUM  3  Weighted  averages  of  the  standard  errors 
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TABLE  5.1.6.  MINIMA  AND  MAXIMA  OF  THE  STANDARD  ERRORS  OF  THE  LOGARITHMIC 
50,  80,  AND  90  PERCENTILES  AND  THE  SLOPE  OVER  ALL  THE 
DISTRIBUTIONS  IN  THE  SENSITIVITY  ANALYSIS 


Minimum  and  maximum  of  unweighted  standard  errors  of  points  varied 
around  A  and  B  over  the  designs  in  the  sensitivity  analysis 

Design  I.D.  3  GDI 


VARNAME 

N 

MINIMUM 

MAXIMUM 

SEFVX50 

49 

0.029316 

0.047555 

SEFVX80 

49 

0.041578 

0.105432 

SEFVX90 

49 

0.055582 

0.144400 

SEFVXB1 

49 

0.805459 

0.919895' 

Design  I.D.  3  GD10 

SEFVX50 

49 

0.031554 

0.047595 

SEFVX80 

49 

0.041651 

0.087573 

SEFVX90 

49 

0.053014 

0.117208 

SEFVXB1 

49 

0.641773 

0.799141 

Design  I.D.  3  GD2 

SEFVX50 

49 

0.030646 

0.047480 

SEFVX80 

49 

0.041790 

0.092691 

SEFVX90 

49 

0.054465 

0.125116 

SEFVXB1 

49 

0.690566 

0.861523 

Design  I.D.  3  GD3 

SEFVX50 

49 

0.031852 

0.047782 

SEFVX80 

49 

0.042300 

0.086001 

SEFVX90 

49 

0.054073 

0.114733 

SEFVXB1 

49 

0.627240 

0.822425 

Design  I.D.  3  GD4 

SEFVX50 

49 

0.030076 

0.047356 

SEFVX80 

49 

0.041038 

0.097405 

SEFVX90 

49 

0.053542 

0.132328 

SEFVX81 

49 

0.733219 

0.849374 
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TABLE  5.1.6. 
(Continued) 


Design 

I.D.  =» 

GD5 

VARNAME 

N 

MINIMUM 

MAXIMUM 

SEFVX50 

49 

0.029622 

0.047440 

SEFVX80 

49 

0.041387 

0.101758 

SEFVX90 

49 

0.054838 

0.138902 

SEFVXB1 

49 

■ 

0.772841 

0.894420 

■ 

Design 

I.D.  * 

GD6 

SEFVX50 

49 

0.02775 

0.04980 

SEFVX80 

49 

0.04441 

0.14178 

SEFVX90 

49 

0.06528 

0.19805 

SEFVXB1 

49 

1.12243 

1.19243 

Design 

I.D.  = 

GD7 

SEFVX50 

49 

0.032382 

0.048034 

SEFVX80 

49 

0.042927 

0.083840 

SEFVX90 

49 

0.054734 

0.111264 

SEFVXB1 

49 

0.605445 

0.826605 

Design 

1.0.  * 

GD8 

SFFVX50 

49 

0.031101 

0.047525 

SEFVX80 

49 

0.041628 

0.089953 

SEFVX90 

49 

0.053529 

0.120888 

SEFVXB1 

49 

0.0664329 

0.824324 

Design 

I.D.  = 

GD9 

SEFVX50 

49 

0.034806 

0.048910 

SEFVX80 

49 

0.044264 

0.077431 

SEFVX90 

49 

0.054499 

0.100846 

SEFVXB1 

49 

0.539006 

0.755271 

VARNAME  Identifier  of  quantity  being  estimated 

N  =  Numbe;’  of  distributions 

entering 

into  the  estimate 

MINIMUM  -■  Minimum  standard 

error  of  the  logarithm 

MAXIMUM  =  Maximum  standard 

error  of  the  logarithm 
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estimate  would  be  updated  and  the  sensitivity  analysis  would  be  carried  out, 
as  above,  to  determine  which  second-stage  dose  allocations  best  augment  the 
first-stage  doses,  in  light  of  what  has  been  learned  about  the  dose-response 
distribution  from  the  first  stage.  This  procedure  would  be  iterated  following 
each  stage  of  experimentation. 

5.2  Probit  Dose  Response  Estimation  Based  on 
Nonlinear  Regression  Analysis 

Following  each  stage  of  experimentation,  the  estimates  of  the 
underlying  dose-response  distributions  are  updated.  Probit  dose-response 
models  in  logarithmic  dose  (Finney  1977)  are  fitted  to  the  data  for  each 
treatment  regimen  to  quantify  the  relationships.  Distribution  percentiles  are 
estimated  based  on  these  models.  Background  response  is  not  incorporated  into 
the  models,  due  to  the  relatively  short  durations  of  the  tests  (hours,  days, 
or  at  most  one  or  two  weeks). 

Standard  probit  analysis  computer  programs  cannot  be  used  to  fit 
these  models  to  the  dose-response  data  due  to  the  nonstandard  dose-allocation 
strategy  and  due  to  a  number  of  nonstandard  aspects  of  the  model 
specifications.  These  nonstandard  aspects  include  individual  animal  responses 
rather  than  pooled  group  lethality  rates,  common  probit  slopes  shared  by 
several  treatment  regimens,  the  possible  presence  of  stage  effects  and  the 
capability  to  adjust  for  such  effects,  and  the  incorporation  of  covariates, 
such  as  body  weight,  into  the  models.  These  model  aspects  are  discussed  in 
greater  detail  in  this  section. 

Specialized  procedures,  based  on  nonlinear  regression  analysis,  have 
been  developed  to  fit  dose-response  models  to  such  data.  These  procedures  are 
described  and  illustrated. 

5.2.1  Individual  Animal  Responses 

The  dose-allocation  strategy  discussed  in  the  previous  section 
results  in  many  different  dcr,es  with  few  animals  tested  per  dose,  possibly 
just  one.  The  model  fitting  methods  thus  need  to  accommodate  the  possibility 
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of  dose-response  data  with  each  animal  tested  at  a  unique  dose.  All  the 
responses  would  then  be  0's  and  l's.  This  is  in  contrast  to  the  usual  probit 
analysis  situation  where  multiple  animals  are  tested  at  a  relatively  small 
number  of  repetitive,  discrete  doses  and  dose-response  models  are  fitted  to 
the  observed  response  proportions  at  each  dose. 

5.2.2  Separate  Slopes  and  Common  Slopes  Models 

The  model  fitting  procedures  are  sometimes  used  to  compare 
dose-response  distributions  corresponding  to  several  treatment  regimens.  For 
example,  no  treatment,  standard  treatment,  and  one  or  more  candidate 
treatments  may  be  compared  simultaneously.  A  fully  general  model  fits 
separate  dose-response  distributions  to  each  regimen.  Submodels  incorporating 
the  assumption  of  common  slopes  among  various  subsets  of  the  treatment  afford 
the  possibility  of  substantially  greater  estimation  precision  and 
interpretation  simplicity.  Provisions  have  been  incorporated  into  the  model 
fitting  procedures  to  fit  common  slopes  to  various  subsets  (of  size  2  to  5)  of 
the  treatments  and  to  test  the  adequacy  of  fit  of  the  submodels  relative  to 
the  separate  slopes  model  or  to  less  restrictive  common  slopes  models.  Common 
slopes  models  are  also  sometimes  used  to  augment  information  about  the  dose 
response  for  a  current  treatment  with  that  based  on  historical  data.  Although 
the  dose-response  distributions  may  be  shifted  relative  to  one  another,  the 
slopes  may  have  remained  the  same. 

Tests  of  adequacy  of  the  submodels  are  carried  out  by  comparing  the 
values  of  the  log  likelihoods  under  the  more  restrictive  and  the  less 
restrictive  models.  The  log  likelihood  tatio  is  referred  to  the  upper 
percentiles  of  a  chi-square  distribution  with  an  appropriate  number  of  degrees 
of  freedom. 

5.2.3  Stage  to  Stage  Variation 

The  basic  design  strategy  calls  for  carrying  out  the  dose-response 
experiment  in  stages,  utilizing  the  results  from  all  previous  stages  to  design 
the  following  stage.  A  test  for  the  presence  of  stage  to  stage  variation  is 
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incorporated  into  the  model  fitting  procedures.  The  test  is  carried  out  by 
fitting  probit  models  to  the  combined  data  across  all  stages.  Residuals  from 
these  fits  are  standardized  by  dividing  them  by  estimates  of  their  standard 
deviations.  In  the  absence  of  stage  to  stage  variation,  these  standardized 
residuals  would  be  expected  to  have  mean  approximately  0  and  standard 
deviation  approximately  1.  If  systematic  stage  to  stage  variation  exists, 
then  the  residuals  from  some  stages  would  have  positive  means  while  those  from 
other  stages  would  have  negative  means. 

A  one-way  analysis  of  variance  is  carried  out  on  the  standardized 
residuals,  incorporating  stage  as  a  grouping  variable  and  the  logarithm  of 
dose  as  a  covariate.  Statistical  significance  of  the  stage  factor  (a  =  0.05) 
provides  evidence  of  stage  to  stage  variation.  Possible  causes  of  such  stage 
effects  might  be  drift  across  stages,  isolated  outlying  responses,  or 
variation  of  some  of  the  experimental  conditions  across  stages. 

The  nature  of  the  stage  to  stage  variation  would  need  to  be  studied 
by  more  in-depth  examination  of  the  data,  such  as  diagnostic  plots,  multiple 
comparison  procedures,  or  the  incorporation  of  additional  explanatory 
variables  into  the  models.  The  nature  and  extent  of  such  additional  analyses, 
and  possible  actions  taken  as  a  result,  would  necessarily  be  decided  upon  on  a 
case  by  case  basis.  They  are  not  incorporated  into  the  more  general  model 
fitting  procedures  discussed  here. 

5.2.4  Covariates 

Body  weight  (kg)  at  the  time  of  dosing  is  incorporated  into  the 
model s  as  a  covariate.  Models  incorporating  separate  covariate  effects  for 
different  treatment  regimens  and  models  incorporating  common  covariate  effects 
are  fitted  to  the  data.  Likelihood  ratio  tests  for  common  covariate  effects 
are  carried  out  in  the  same  manner  as  likelihood  ratio  tests  for  common 
dose-response  slopes. 
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It  should  be 
the  dose-response  model 
responses  rather  than  p 
dose.  This  is  because, 
and  so  presents  a  diffeh 


5.2.5  Dose-Response  Model  Fitting  Procedures 


rioted  that  including  covariates  such  as  body  weight  in 
necessitates  fitting  to  the  individual  animal  0-1 
<j>oling  across  animals  that  were  tested  at  the  same 
in  general,  each  animal  has  a  different  body  weight 
ent  set  of  explanatory  variables. 


A  series  of  computer  programs,  based  on  PROC  NUN  in  the  SAS 
statistical  computing  system,  have  been  developed  to  fit  the  dose-response 
models  to  the  experimental  results.  These  procedures  utilize  as  input  the 
individual  animal  0-1  responses,  as  well  as  the  treatment  dose  and  any 
covariates,  such  as  body  weight.  Programs  are  available  to  fit  separate 

I 

probit  models  to  each  individual  treatment  (separate  slopes  model)  and  to  fit 
joint  probit  models  having  a  common  slope  to  several  treatments  (common  slopes 
model).  Covariates  can  be  included  in  or  excluded  from  the  models. 

Table  5.2.1  displays  the  output  from  a  common  slopes  probit  model 


fit  to  the  results  from! 


two  treatment  regimens.  The  parameter  81  represents 


the  common  slope  and  the  parameters  801,  802  represent  the  intercepts  for  the 
treatments.  No  covariaite  is  included  in  this  model.  If  the  model  fits  the 


data,  then  the  expected 


value  of  the  residual  mean  square  is  asymptotically 


1.0.  The  attained  residual  mean  square  of  0.85  indicates  no  evidence  of  lack 
of  fit  of  the  model.  The  "sum  of  loss"  is  proportional  to  -2  times  the 
(natural)  logarithm  of  [the  likelihood  function;  it  is  used  to  compare  the 
adequacy  of  alternative  models. 

Tables  5.2.2  and  5.2.3  display  estimates  of  the  dose-response 
distribution  percenti lejs,  associated  standard  errors,  and  upper  and  lower 
95  percent  confidence  bounds.  The  confidence  bounds  in  Table  5.2.2  are  based 
onpropagation  of  errors,  while  those  in  Table  5.2.3  are  based  on  Fieller's 
method.  If  the  estimated  slope  is  somewhat  more  than  two  standard  errors 
from  0,  as  in  this  example,  then  both  confidence  intervals  are  similar, 
particularly  for  doses  iin  the  central  portion  of  the  design.  If  the  estimated 
slope  is  less  than  two  standard  errors  from  0,  then  the  Fieller's  method 
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confidence  intervals  will  be  substantially  wider  than  the  propagation  of 
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TABLE  5.2.1.  COMMON  SLOPES  PROBIT  DOSE-RESPONSE  MODEL  FIT  TO  THE  RESULTS 
FROM  TWO  TREATMENT  REGIMENS 


ABLE  5.2.2.  ESTIMATED  DOSE-RESPONSE  DISTRIBUTION  PERCENTILES  AND  ASSOCIATED 
CONFIDENCE  BOUNDS  BASED  ON  PROPAGATION  OF  ERRORS 
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TABLE  5.2.3.  ESTIMATED  DOSE -RESPONSE  DISTRIBUTION  PERCENTILES  AND  ASSOCIATED 
CONFIDENCE  BOUNDS  BASED  ON  FIELLER'S  HETHOD 
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errors  intervals;  the  propagation  of  errors  intervals  are  too  narrow,  while 
the  Fieller's  method  intervals  are  too  wide.  A  compromise  interval  cannot  be 
obtained  analytically;  it  likely  requires  a  resampling  method,  such  as 
bootstrapping,  to  account  for  the  inherent  nonlinearities. 

Tables  5.2.4  and  5.2.5  display  summaries,  by  stage,  of  the 
standardized  residuals  from  the  probit  model  fit  and  a  one-way  analysis  of 
variance  to  test  for  the  presence  of  stage  to  stage  variation  in  these 
residuals.  There  is  no  evidence  of  a  significant  stage  effect  in  this 
example.  If  the  stage  effect  in  Table  5.2.5  was  significant,  then  the 
stagewise  means,  standard  deviations,  minima,  and  maxima  in  Table  5.2.4  would 
be  studied  to  determine  the  nature  of  the  variation,  and  which  stage  or  stages 
differ  from  the  remainder. 

Tables  5.2.6  and  5.2.7  display  the  outputs  from  two  probit  model 
fits  to  the  results  from  a  different  dose-response  experiment,  with  two 
treatments.  Body  weight  (kg)  is  included  as  a  covariate  in  these  models. 

Table  5.2.6  displays  the  results  of  a  four-parameter  common  slopes  model,  with 
a  common  covariate  effect,  fitted  to  the  two  treatments.  The  parameters  B1 
and  P2  represent  the  common  slope  and  the  common  body  weight  effects, 
respectively;  B01  and  B02  represent  the  intercepts  corresponding  to  treatments 
1  and  2,  respectively.  Table  5.2.7  displays  the  results  of  a  single  three- 
parameter  probit  model  fitted  to  the  combined  results  from  both  treatments. 

B!  and  82  represent  the  slope  and  the  body  weight  effect,  respectively; 

60  represents  the  intercept.  Based  on  the  residual  mean  square,  both  models 
appear  to  fit  the  data. 

A  log  likelihood  ratio  test  for  differences  between  the  dose- 
response  distributions  is  carried  out  by  comparing  the  difference  between  the 
"sum  of  loss"  values  for  the  two  models  to  a  chi-square  distribution  with 
1  degree  of  freedom  (4  parameters  minus  3  parameters).  Namely,  60.989  - 
60.538  *  0.451  is  significant  at  the  a  *  0.50  level,  based  on  the  chi-square 
distribution  with  1  degree  of  freedom.  Thus,  there  is  no  evidence  of 
differences  between  the  dose-response  distributions  associated  with  each  of 
the  treatments  and  so  the  single  model  is  accepted. 


I  Ad.,  b  2.4.  iuatAKitS.  BY  STAGE.  Of  IH£  STANDARDIZED  RESIDUALS  FROM  THE  PR061T  MODEL  FIT 
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The  body  weight  parameter  can  be  interpreted  as  follows.  The 
logarithm  of  the  LDSI  dose  for  a  W  kg  animal,  based  on  the  model  in 
Table  5.2.7,  can  be  calculated  by  solving  the  equation 

BO  -  5  +  B1  *  x5i  +  B2(W  -  W)  *  0, 

where  xsl  represents  the  estimated  common  logarithm  of  the  LDSi  dose  and  W  is 
the  average  body  weight.  Thus, 

A  BO-5  B2 

x5>  3  -  - -  -  -  (W  -  W) . 

B1  81 

If  B2  is  positive,  then  this  relation  can  be  interpreted  as  a  decrease  in  the 
LDSI  of  100  (1-10  '82/B1)  percent  tor  each  1-kg  increase  in  body  weight.  For 
the  present  example,  this  is  an  estimated  2.6  percent  decrease  in  the  LDSI  for 
each  kg  increase  in  body  weight,  at  least  for  body  weights  around  the  average. 

6.0  POSSIBLE  DIRECTIONS  FOR  FURTHER  EXTENSIONS 

This  section  considers  several  possible  directions  for  further 
development,  extension,  or  modification  of  the  methods  ana  procedures  that 
were  discussed  in  the  previous  sections. 

6.1  Parameter  Selection 


The  primary  reasons  for  the  selection  of  the  input  parameter  values 
of  the  recommended  test  procedure  are  that  these  values  are  consistent  with 
the  historical  database  and  they  appear  to  provide  a  test  procedure  with 
desirable  overall  properties  as  illustrated  in  Sectipn  2.3.  Further  work  must 
be  performed  to  develop  a  procedure  for  the  selection  of  test  procedure 
parameters,  allowing  these  parameters  to  vary  with  the  agent  and  test  system. 
The  two  major  activities  that  would  be  required  are  more  extensive  simulation 
studies  to  characterize  the  behavior  of  the  recommended  test  procedure  for 
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various  input  parameter  values  and  a  statistical  analysis  of  the  historical 
database,  using  the  methods  of  Section  3.2,  to  determine  plausible  values  of 
the  input  parameters. 

6.2  Historical  Data 

The  methods  discussed  in  Sections  2.0  and  3.0  assume  that  an 
extensive  amount  of  historical  data  is  available.  The  parameters  nc,  o^, 
and  e  that  characterize  the  historical  response  distribution  for  the  standard 
decontaminant  are  assumed  known.  In  some  applications,  however,  there  might 
be  just  a  small  or  moderate  amount  of  historical  data  available.  The 
historical  distribution  parameters  would  then  have  uncertainty  associated  with 
their  estimates.  This  uncertainty  would  inflate  the  variability  of  the 
historical  estimates  relative  to  the  expressions  presented  in  Section  2.0  and 
would  thereby  result  in  greater  weight  being  given  to  the  current  estimate. 

The  methods  discussed  in  this  report  can  be  extended  to  account  for  this  ■ 
additional  source  of  variability  and  its  influence  on  the  recommended 
weighting  procedure. 


6.3  Deleting  Far  Past  Historical  Data 

Procedures  might  be  developed  for  determining  when  and  to  what 
extent  to  delete  the  far  past  standard  decontaminant  results  when  they  are  no 
longer  compatible  with  the  current  and  more  recent  past  standard  decontaminant 
results.  Such  decisions  would  be  based  on  exceedences  observed  with  the 
control  chart  procedures  discussed  in  Section  3.0. 

6.4  Discounting  Historical  Data  Based  on  Its  Age 

Current  procedures  utilize  all  the  historical  data  as  equivalent  as 
long  as  they  remain  within  the  control  limits  and  as  long  as  weighted  averages 
of  various  durations  remain  within  the  control  limits.  This  is  the  case 
whether  the  historical  values  were  obtained  a  day,  a  week,  or  a  year  ago.  An 
alternative  procedure  is  to  routinely  discount  the  historical  data  based  on 


their  age,  irrespective  of  whether  or  not  they  lie  within  the  control  limits. 


For  example,  if  exponential  discounting  were  used  w 
-X  per  month,  then  one  current  observation  would  be 
observations  in  a  month,  e'12^  observations  in  a  yearj 
screen  were  not  used  on  a  regular  basis,  the  extent 
information  would  gradually  diminish  over  time.  If 

historical  information  about  the  standard  decontamiijiant  response  rate  drops 
below  a  specified  level,  then  additional  tests  would  be  carried  out  with  the 
standard  decontaminant  to  increase  the  amount  of  historical  information  up  to 
a  specified  minimum  threshold.  These  additional  tests  would  be  carried  out  at 
the  same  agent  dose  as  that  used  in  past  tests.  If  the  standard  decontaminant 
response  rates  drift  (or  jump)  out  of  control,  as  determined  by  the  control 


th  a  discount  parameter  of 
discounted  to  e*^ 

,  etc.  Thus,  if  the 
of  the  historical 
the  amount  of  discounted 


chart  procedures,  then  a  new  LDSI  study  would  be  car 
to  modify  the  agent  dose  for  future  tests. 


6.5  Beta  Binomial  Distribution 


ried  out  to  determine  how 


The  observed  response  rates,  for  both  the  standard  and  test 
decontaminants,  are  currently  modelled  as  being  approximately  normally 
distributed.  Following  an  arc  sin  transformation,  the  variances  of  these 
response  rates  are  assumed  to  be  independent  of  the  mean.  The  normal  , 
approximation  to  the  binomial  distribution  is  reasonable  for  response  rates 
near  50  percent,  as  is  the  case  with  the  current  application.  For  other 
applications,  with  response  rates  closer  to  0  or  100  percent,  the  normality 
assumption  may  not  be  as  appropriate.  An  alternative  formulation  for  such 
problems  would  be  to  model  the  responses  within  each  individual  test  as 
binomial ly  distributed  with  response  probability  vairying  among  tests  according 
to  a  beta  distribution.  The  resulting  marginal  distribution  of  the  observed 
standard  decontaminant  response  rates  across  tests  can  be  described  by  the 
beta  binomial  distribution.  This  distribution  is  bounded  between  0  and  1  and 
incorporates  skewness  in  the  appropriate  direction  when  the  true  response 
probabilities  are  near  0  or  1. 


There  are  no  conceptual  differences  in  a  model  formulation  based  on 
the  beta  binomial  distribution  from  one  based  on  the  normal  distribution. 

There  are,  however,  a  number  of  technical  differences;  the  expressions  for  the 
weights  would  need  to  be  modified. 

6.6  Determination  of  Control  Chart  Boundaries 


In  Section  2.0,  the  test  to  test  variation  of  standard  decontaminant 
response  rates  is  modelled  as  a  mixture  of  normal  distributions.  The  control 
chart  limits  in  Section  3.0,  however,  are  based  on  a  single  normal 
distribution.  The  distributional  assumptions  made  in  Section  3.0,  and  control 
limits  based  on  them,  might  be  modified  to  be  brought  into  conformance  with 
the  assumptions  made  in  Section  2.0.  Namely,  the  control  limits  might  be 
based  on  the  upper  percentiles  of  the  mixture  distribution,  using  the 
iterative  calculation  recommended  in  Section  2.2  to  determine  the  critical 
value.  This  would  probably  result  in  wider  control  limits  than  those  based  on 
the  normal  approximation  to  this  distribution. 

6.7  Determination  of  Control  Chart  Statistics 

The  discussion  in  Section  3.3  refers  to  three  alternative  statistics 
to  indicate  when  the  standard  decontaminant  response  rates  are  drifting  away 
from  historical  levels.  An  individual,  standardized  transformed  response  rate 
is  associated  with  each  test.  One  statistic  is  based  on  the  numbers  of 
consecutive  individual  values  that  exceed  control  limits.  A  second  statistic 
is  based  on  comparing  the  medians  of  rjnsecutive  individual  values  to  control 
limits.  A  third  statistic  is  based  jn  comparing  the  means  of  consecutive 
individual  values  to  control  limits. 

Table  3.4.1  in  Sect'sn  3.4  shows  that  the  statistic  based  on  the 
means  is  more  powerful  for  detecting  small  to  moderate  departures  than  those 
based  on  counts  or  medians.  However,  the  statistic  based  on  the  means  is  more 
sensitive  to  the  effects  of  a  small  number  of  outlying  values.  A  compromise 
between  the  means-based  statistic  and  the  medians-based  statistic  might  be 
found  that  simultaneously  provides  much  of  the  improved  sensitivity  to  detect 


systematic  departures,  yet  resists  much  of  the  insensitivity  and  is  not  as 
influenced  by  isolated  outlying  values.  Such  a  compromise  procedure  might  be 
based  on  trimmed  means  of  consecutive  values. 

6.8  Generalization  of  the  Dose-Response  Models 

The  discussion  in  Sections  4.0  and  5.0  pertains  to  experimental 
design  and  data  analysis  considerations  in  deterrining  the  LDsg  associated 
with  the  standard  decontaminant.  It  is  assumed  there  that  the  dose-response 
relation  can  be  described  by  a  probit  model,  without  background.  This  model 
is  adequate  for  many  applications  to  which  the  screening  methodology  has  been 
applied.  Other  applications,  however,  might  necessitate  the  use  of  more 
general  models. 

A  test  period  of  relatively  long  duration  might  result  in  a  nonzero 
background  lethality  response  rate.  Treatment  with  a  specified  drug  regimen 
might  not  be  efficacious  for  all  the  animals,  no  matter  how  much  the  drug  dose 
is  increased.  Morbidity  responses,  such  as  deterioration  of  neurological 
function,  may  be  exhibited  by  some  animals,  no  matter  how  high  the  drug  dose, 
and  may  not  be  exhibited  by  some  animals,  even  in  the  absence  of  drug 
treatment. 

The  sensitivity  analysis  procedures  and  the  dose-response  model 
fitting  procedures  discussed  in  Sections  4.0  and  5.0  might  be  extended  to 
accommodate  minimum  and  maximum  response  rates  strictly  between  0  and  1. 

These  rates  would  be  additional  model  parameters,  to  be  estimated  from  the 
data. 
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APPENDIX  A 


DOCUMENTATION  FOR  NEW  SAS  PROGRAMS 


DOCUMENTATION  FOR  NEW  SAS  PROGRAMS 
A.l  DATA  STRUCTURES 
A. 1.1  Historical  Data  File 

The  historical  data  file  is  an  ASCII  file  containing  standard 
decontaminant  results  from  first  stage  screens.  There  is  a  separate  file  for 
each  agent.  Each  record  contains  information  about  a  screen,  including  the 
starting  date,  the  number  of  animals  dosed,  and  the  number  of  lethalities. 

The  first  field  in  each  record  is  a  “USEFLAG"  which  can  be  set  to  either  use 
or  ignore  that  record  when  computing  historical  estimates  and  creating  control 
charts.  As  ner  standard  decontaminant  data  becomes  available,  it  can  be 
appended  to  the  end  of  the  existing  historical  data  file. 

File  name:  <agent>.HIS,  where  agent  is  the  agent  code. 


Field 

Name 

Columns 

Description 

USEFLAG 

1 

Flag  which  is  set  to  0  or  1  to  ignore  or  use  data 
when  calculating  the  historical  estimate. 

TSEQ 

3  -  4 

First  stage  screen  identifier  (A,  B,  C...) 

AGTCD 

6  -  8 

Agent  code  (GO,  TGD,  VX) 

DCNCD 

10-12 

Decontamination  code  (STD  *  standard,  A,B,C. . .*test) 

STRTDATE 

14-21 

First  date  of  testing  for  the  screen  (mm/dd/yy 
format) 

DURATION 

24-30 

Duration  of  screen  in  days:  (last  date)  -  (first 
date)  +1 

NOOSED 

32-38 

Number  of  animals  which  were  dosed 

NDEAD 

40-46 

Number  of  animals  which  died 

A-l 


A  1.2  Current  Data  Files 


The  current  data  file  is  an  ASCII  file  containing  results  of  a 
single  first-stage  screen  using  the  standard  system  and  the  test 
decontaminants  for  a  particular  agent.  This  file  is  in  the  same  format  as  the 
historical  data  file  described  above.  Each  record  corresponds  to  a  single 
standard  or  test  :iecontaminant. 

File  name:  <agent>.CUR  where  <agent>  is  the  agent  code. 

A. 1.3  Nominal  Parameter  Values  Files 

One  file  per  agent,  containing  the  nominal  values  for  the  lethality 
rate  and  the  screen-to-screen  variability. 


File  name:  <agent>.N0M  and  <agent>.EST  where  agent  is  the  agent  code. 


Field 

Name 

Columns 

Description 

MUC 

1-20 

Nominal  lethality  rate,  pc. 

SIGD2 

21-40 

Nominal  screen-to-screen  variance,  a 

A. 2 

DOCUMENTATION  FOR  SCREEN  PROGRAMS 

The  following  programs  perform  the  analyses  in  the  attached  Report 
on  Lethality  Rate  Estimation  and  Testing  Procedures.  The  programs  are  written 
using  in  the  Statistical  Analysis  System  (S-‘c)  and  are  designed  to  run  on 
Battelle's  VAX  system.  It  is  assumed  that  c  'ogram  and  data  files  reside 
in  the  default  directory  of  the  analyst. 

To  run  the  programs,  log  on  to  the  Battelle  computing  network 
requesting  a  destination  of  VMSF,  and  at  the  VMS  prompt  (F$),  type: 

SAS  <program  name> 


A-2 


followed  by  a  carriage  return.  Each  program  produces  a  file  named  <program 
name>.LIS  which  contains  the  analysis  results,  and  a  file  named  <p^ogram 
name>.L0G  which  contains  the  SAS  log. 

A. 2.1  Documentation  for  "HISTLETH" 

Program  file:  HISTLETH. SAS 

This  program  uses  a  supplementary  file  named  HISLETHS. SAS 

This  program  computes  estimates  of  lethality  rate  and  screen-to- 
screen  variance  using  historical  data  from  first-stage  screens.  The  required 
input  to  the  program  is  a  file  called  <agent>.HIS  which  contains  historical 
data  on  the  standard  system  for  a  particular  agent.  Output  of  the  program 
consists  of  a  printed  report  and  a  data  file.  The  printed  report  includes  the 
input  data,  estimates  of  lethality  rate  and  screen-to-screen  variability,  and 
standard  errors  for  these  estimates.  The  data  file  is  called  <agent>.EST 
and  includes  the  estimates  of  lethality  rate  and  screen-to-screen  variance. 

The  format  of  this  file  is  given  in  A. 1.3. 

Prior  to  running  the  program,  edit  HISTLETH. SAS  to  include  the  names 
of  the  .HIS  historical  data  file  to  be  used  and  the  . EST  file  to  be  created. 

A  listing  of  the  HISTLETH  program  is  provided  in  Section  8-1  of 
Appendix  B.  A  listing  of  the  HISLETHS  file  is  provided  in  Section  8-2. 

A.2.2  Documentation  for  “COMPARE" 


Program  file:  COMPARE. SAS 

This  program  uses  supplementary  files  named  CQNTAMC.SAS  and 
CONTAMS . SAS .  _  . 

This  program  carries  out  the  statistical  test  procedures  for 
comparing  the  current  behavior  of  the  standard  decontaminant  with  historical 
behavior,  and  for  comparing  each  experimental  decontaminant  with  the  standard 
•iecontominant.  Two  data  files  are  required  as  input  to  the  program.  The 
first  is  a  tile  called  •'agents .CUR  which  contains  current  test  system  data 
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including  standard  system  and  test  decontaminant  results.  The  second  is 
called  <agent>.N0M  and  contains  estimates  of  the  nominal  historical  lethality 
rate  and  screen-to-screen  variance.  Output  of  the  program  consists  of  a 
report  and  a  data  file.  The  printed  report  contains  listings  of  the  raw  data 
and  computed  values  and  significance  levels  from  the  test  procedures.  The 
data  file  is  named  <agent>.NEWHIS  and  is  an  ASCII  file  containing  the  current 
standard  decontamination  data.  This  file  can  be  APPENDed  to  the  historical 
data  file  by  the  analyst. 

Prior  to  running  the  program,  edit  COMPARE. SAS  to  include  the  names 
of  the  .CUR  current  system  data  file  and  the  .NOM  file  to  be  used,  and  the 
.NEWHIS  data  file  to  be  created. 

A  listing  of  the  COMPARE  program  is  provided  in  Section  B-3  of 
Appendix  B.  Listings  of  the  CONTAMC  and  CONTAMS  files  are  provided  in 
Sections  B-4  and  B-5,  respectively. 

A. 2. 3 _ Documentation  for  "CRITX" 


Program  file:  CRITX. SAS 

This  program  uses  supplementary  files  named  CRITX.DAT  (which 
must  be  created  by  the  user),  CONTAMC. SAS  and  CONTAMS. SAS. 

This  program  determines  the  critical  values  for  the  test  procedure 
comparing  an  experimental  decontaminant  with  the  standard  decontaminant.  The 
analyst  specifies  values  of  certain  parameters  and  the  program  produces  a 
table  of  critical  values  for  determining  when  a  test  decontaminant  is  no' 
better  than  the  standard  decontaminant  based  on  the  number  of  lethalities  in 
the  standard  and  test  groups.  The  six  parameters  specified  by  the  analyst 
are: 

nc  number  of  animals  receiving  the  standard  decontaminant 

nt  number  of  animals  receiving  the  test  decontaminant 

Hc  long-term  lethality  rate  for  the  standard  decontaminant 

e  mixture  probability  for  the  random  effect  5 
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<7 standard  deviation  of  nominal  normal  distribution 
standard  deviation  of  extreme  normal  distribution 

Values  for  these  six  parameters  should  be  entered  into  a  one  line  ASCII  file 
named  CRITX.DAT  in  the  order  specified  above,  with  at  least  one  blank  space 
separating  each  value. 

The  program  produces  a  file  named  CRITX.LIS  displaying  the  input 
parameters,  and  a  table  of  critical  values.  The  critical  values  are  tabulated 
for  values  of  xc,  the  number  of  lethalities  for  the  standard  decontaminant, 
ranging  from  0  to  nc.  This  information  is  also  printed  on  the  screen  as  the 
program  runs. 

A  listing  of  the  CRITX  program  is  provided  in  Section  B-6  of 
Appendix  B.  Listings  of  the  CONTAMC  and  CONTAMS  files  are  provided  in 
Sections  B-4  and  8-5 ,  respectively. 
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8-1  PROGRAM  LISTING  FOR  "HISTLETH 


DATA  HISTORIC; 


SUPPLY  THt  NAM  OP  THC  HISTORICAL  OATA  PXLf  IN  TH«  POLLOVINQ 
INPILE  STATEMENT 


INPILE 

'SO. HIS' ; 

INPUT 

USEPLAQ 

1 

n  tseq 

$2. 

99  AST CO 

S3. 

P10  OCNCO 

S3. 

•14  STRTOATE 
P24  DURATION 
*22  NOOSED 
*40  NOCAO: 

MNOOYYS . 

PORMAt  STRTOATI  MNOOYYS.: 

PROC  PRINT;  TITLE  'HISTORICAL  STANOARO  OATA' : 
• 

•  GENERATES  HISTORICAL  ESTIMATES 

a  • 

OATA  OUMIV; 

SET  HISTORIC; 

IP(USEPLAO  EQ  1); 

OPTIONS  NOSOURCE2; 

XINCLUOE  HISLETHS: 

OATA  HISTEST; 

SET  CURREST; 


.  "+* 


SUPPLY  THE  NAME  OP  THE  NOMINAL  VALUE  PILE  IN  THt  POLLQWINQ 
PILE  STATEMENT 


PILE  'QO. EST' ; 

PUT  *1  PH  *21  SI002H: 

PROC  PRINT;  TITLE  'HISTORICAL  ESTIMATES'; 
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B-2  PROGRAH  LISTING  FOR  “HISLETHS" 


OATA  CURREST; 

INPUT  PH  SI GOSH; 

CAROS; 

0.9  0.0 

i  ARCSIN-SORT  TRANSFORMATION  is  USED  TO  TRANSFORM  THt  BINOMIAL 
.  NOOIL  to  constant  VAAIANCI 


•  ITERATION  1 

•  ; 

OATA  ITER; 

SET  ourar ; 

IF  N  EQ  1  then  SET  CURREST; 

Wd*7/7. 2S/N0QSED+SIG02H) ; 

Zd* 1 / ( 2»  < • 2S/NOOSEO*SIG02H) *  *2 ) ; 

P o  •  ARSIN(SORT(NOEAO/NOOSEO) )  ; 

SIQ2J  *(Pd-PH)**2  -  . 23/NOOSED; 

HUMERI  •  Wd»Pd: 

NUMER2  ■  Zd»SI02d; 

PROC  SUMMARY  OATA* ITER; 

VAR  NUMERI  NUMER2  Wd  Zd; 

OUTPUT  OUT* SUM 1  5UM*SUMPd  SUMSIQ2J  SUMWd  SUMZd; 

OATA  CURREST; 

SET  CURREST; 

IP  <_N_  EQ  1)  THEN  SET  SUM1 ; 

PP«*pn7 

pa 1902ft* a 1902ft: 

PH  *  SUMPO/ SUMWd; 

SEPH  •  SQRTt 1 /SUMWd): 

SI002H*SUMSIQ2d/SUMZd ; 

SESIG02H*$QRT( 1/SUMZd); 

PCPH  •  100* ( PH-PPH ) /PPH ; 

PCSI002H  •  100* ( SIG02H-PSIG02H ) /PSIG02H; 

PUT  ' PCT  CHANGE  IN  PH  »  '  PCPH  '.  PCT  CHANGE  IN  SIQ02H  »  ' 
KEEP  PH  SEPH  SIG02H  SESI002H; 

« 

•  ITERATION  2 
» 

OATA  ITER: 

set  ouirr. 

IF  N  EQ  1  THEN  SET  CURREST; 

Wd*7/7. 2S/N00SE0*SIG02H) ; 

Zd*1/( 2*( . 23/N00SED*SIG02H)**2  ) ; 

Pd  ■  ARSINI SQRT(NOEAO/NOQSEO) )  ; 

SI32d  »(Pd-PH)**2  -  . 2S/NOOSEO; 

HUMERI  *  Wd*Pd; 

HUMER2  *  Zd«SIQ2d ; 

PROC  SUMMARY  DATA* ITER; 

VAR  HUMERI  NUMER2  Wd  Zd: 

OUTPUT  OUT* SUM 1  SUM*SUMPd  SUMS I 02 d  SUPWd  SUMZd: 

OATA  CURREST; 

SET  CURREST; 

IF  (  N  EQ  1)  THEN  SET  SUM1 ; 
ppft*pr«7 

palgd2ft*« lg02ft: 

PH  •  SUMPd/ SUMWd; 

SEPH  •  SORT (1/ SUMWd); 

SI002H*SUMSIQ2d/ SUMZd ; 

SESI002H*SQRT( 1/SUMZd ) ; 

PCPH  •  100* (PH-PPH) /PPH; 

PCSI002H  *  100*(SIC02H-PSI002H)/PSI002H; 

PUT  'PCT  CHANQE  IN  PH  •  '  PCPH  '.  PCT  CHANGE  IN  SI002H  •  ' 
KEEP  PH  SEPH  SIS02H  SESIG02H; 

• 

•  ITERATION  3 

a 

OATA  ITER: 

SET  DUMMY ; 

IF  N  Eq’i  THEN  SET  CURREST; 

Wd»7/7. 2S/NOOSEO*SIG02H ) ; 

ZJ* 1/ ( 2*(  . 23/NQOSE0+SIGO2H ) • ■ 2 ) ; 

Pd  *  ARSINI SORTf HOSAO/NOOSEO))  : 

SI02d  *(Pd-PH)«*2  -  .  2S/N00SE0: 

HUMERI  •  Wd*Pd; 

HJMER2  *  Zd»SIG2d: 

PROC  SUMMARY  OATA* ITER; 

VAR  HUMERI  NUMER2  Wd  ZJ; 


PCSIG02H; 


PCSIG02H; 
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OUTPUT  OUT* SUM 1  SUM*SUMPU  SUMS  1 02  d  SUMWd  SUMZJ; 

OAT*  CURREST ; 

SIT  CURRIST; 

IP  <_N_  IQ  1)  THIN  SIT  SUM1; 

PPH*P"7 

pRigo2h*«igaan: 

PH  *  SUMPU/ SUMWd ; 

SIPH  •  SQAT  (  1  /  SUMWvJ ) : 

SI002H*SUMSia2d/SUMXd ; 

SFSIC02H*SQRT( 1/SUMZd): 

PCPM  •  100* ( PH-PPH (/PPH: 

PCSI002H  •  100* ( 3IG02H-PSIQ02H ) /PSIQ02H: 

PUT  '  PCT  CHANOI  IN  PH  •  '  PCPH  '  .  PCT  CHANOI  IN  SIQOSH  • 
KIIP  PH  SIPH  SIQ02H  SISI002H; 


•  ITERATION  4 

• 

OATA  ITER: 

SET  QUMMT 

IP  N  EQ  1  THEN  SIT  CURREST; 

Wd*7/T.  2S/NOOSE0«SISO2H ) : 

Zd*l/(2*( . 2S/N00SID*SIG02H ) • *2 ) : 

Pd  *  *RSIN(SQRT(NOEAO/NOOSEOn  : 

SI02d  »(Pd-PH)**S  -  . 2S/N00SED; 

NUMERI  •  Md»Pd ;  I 

numers  •  Zd«siaaJ; 

PROC  SUMMARY  OATA* ITER; 

VAR  NUMIRI  NUMIRS  Wd  2d; 

OUTPUT  OUT*  SUM  1  SUM*SUMPd  SUMS  1 02  d  SUMWd  SLIMZd; 

OATA  CURREST: 

SIT  CURREST; 

IP  (  N  EQ  1)  THIN  SET  SUM1 ; 

ppn*p«T 

P«iga2n*«igd2h; 

PH  ■  SUMPd/SUMWdJ 
SIPH  •  SORT! 1 /SUMWd ) : 

SI002H*SUM3I02d/SUMZd; 

SESI002H*SQRT( 1/$UMZd> ; 

PCPH  •  100*  <  PH-PPH ) / PPH : 

PCS I S02H  *  100*< SI002M-PSIG02H)/PSIG02M: 

PUT  PCT  CHANOI  IN  PH  •  '  PCPH  '  PCT  CHANOI  IN  SI802N  • 
KEEP  PH  SEPH  SI002H  SESIG02H; 


PCSI002H; 


PCSIQ02H; 


•  ITERATION  S 


OATA  ITER: 

SIT  DU  WIT . 

IP  N  EO  1  THEN  SET  CURREST: 

Wd*1/7.  M/NOOSEO-fSIGOSH) ; 

2d*1/(2*(  .23/N00SE0-*Sia02H)«*2) ; 

Pd  «  ARSIN( SQRT(NOEAO/NOOSEO) )  ; 

SIQSd  *<Pd-PH)«*2  -  . S9/NOQSE0; 

NUMERI  •  Wd*Pd; 

NUMIRS  •  Id*SIQ2d ; 

PROC  SUMMARY  OATA* ITER; 

VAR  NUMERI  NUMERI  Vd  Zd; 

OUTPUT  OUT* SUM 1  $UM*SUNPd  IJMSIQSd  SUMWd  SUMZd: 

OATA  CURREST; 

SET  CURREST: 

IP  (_N_  EQ  1)  THEN  SIT  SUM1 ; 

pph*ph7 

pstgd2n*«1gd2h; 

PH  *  SUMPd/ SUMWd: 

SEPH  *  SQRT( 1/SUMWd ) ; 

SI00SH*SUMSI02d/SUMZd ; 

SE3IC02H*SQRT( 1/SUMZd) ; 

PCPH  •  100* ( PH-PPH ) /PPH; 

PCSICOSH  *  100*(SI002H-PSI002H)/PSI502H: 

PUT  PCT  CHANOI  IN  PH  ■  '  PCPH  PCT  CHANOI  IN  SIG02H  » 
KEEP  PH  SEPH  SIG02H  SESIS02H; 


PCSIS02H: 


*  ITERATION  9 

m 

DATA  ITER: 

SET  DUMMY; 

IP  N  EQ  1  THEN  SET  CURREST; 
Wd* 1/T. 2S/N003E0+SI002H ) ; 


Zd*  1  / ( a* <  . 2S/N00SED*SIG02H) «*J  ) ; 

PU  ■  ARSIN( SQRT( NCEAO/NOOSED ) )  ; 

SI03J  *<Pd-PMI*»2  -  .  25/ NOOSED; 

NONE A 1  •  Wd«Pd: 

NUMER2  •  Zd*SXG2d; 

PPOC  SUMMARY  DATA* ITER ; 

VU  NUMER1  NUMER2  Wd  Zd: 

OUTPUT  OUT* SUM 1  SUM*SUMPd  SUMS I 02 d  SUMWd  SUMZd; 

DATA  CUPP 1ST; 

SIT  CURPIST; 

IP  (  N  IQ  1)  THIN  SIT  SUM1 ; 

PP"«P*7 

P«t9d2ft*ii9d2h: 

PH  «  SUNPd/SUMWd: 

SEPH  •  SORT ( 1 /SUMWd ) : 

SXG02H*SUMSXQ2d/ SUMZd; 

SISI002H*SQRT( 1/SUMZd ) ; 

PCPH  •  100* <  PH-PPHI/PPH; 

PCSXQ02H  •  100*<  SI002H-PSXG02H)/PSXGD2M; 

PUT  'PCT  CHANGE  IN  PM  •  '  PCPH  PCT  WANGE  IN  SIG02H  • 
KEEP  PH  SEPH  SIG02H  SESXG02H; 


PCS I GO 2H; 


•  ITERATION  7 

* 

DATA  ITER; 

set  ouwrr 

IP  N  EQ  1  then  SET  CURREST; 

Wd*?/7.2»/N00SI0-*SIG02H>; 

ZJ»1/( 2*( . 2S/NOOSE0*SIGO2H) ««2 ) ; 

Pd  •  ARSXN(SQRT(NOEAO/NOOS£D) )  ; 

SI62d  »(Pd-PH)«*2  -  . 23/NOOSEC; 

NUMER1  *  Wd*Pd: 

NUMER2  •  Zd*SIS2d ; 

P«OC  SUMMARY  OATA*ITER; 

VAR  NUMSRI  NUMIR2  Wd  Zd; 

OUTPUT  OUT* SUM 1  SUM*SUMPd  SUMS I 02 d  SUMWd  SUMZd; 

DATA  CURREST; 

SIT  CURREST; 

IP  <_N_  EO  1)  THEN  SIT  SUM1 ; 

PP"*P«7 

pstgd2H*f lgM2h; 

PH  •  SUMPd/ SUMWd; 

SEPH  «  SQRT( 1 /SUMWd ) ; 

SIG02H*SUMSIQ2d/ SUMZd; 

SE$IG02H*SQRT( 1 /SUMZd I ; 

PCPH  «  100* (PH-PPHI/PPH: 

PCS I 002 H  *  100*<SIG02H-PSIG02H)/PSIG02H; 

PUT  'PCT  CHANGE  IN  PH  «  '  PCPH  PCT  CHANGE  IN  SX002H  • 
KEEP  PH  SEPH  SXQ02H  SESXS02H; 

*  TRANSPORN  RACK  TO  ORIGINAL  UNITS 
■ 

DATA  CURREST: 

SET  CURREST; 

SEPH  «  2*3IN( PH) *COS( PH) » SEPH; 

PH  •  SIN<  PH)*«2; 


PCS I GO 2H: 
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8-3  PROGRAM  LISTING  FOR  "COMPARE 


DATA  HIST; 


SUPPLY  THt  NAM 


OF  THt  NOMINAL  VAUJt  FILS  IN  THt  FOLLOW I NO 
INF Z LI  STATIMINT 


INFILI  'QO. NOM' ; 
INFLIT  MUC  SI 002  ; 
OATA  STOOATA  Tt STOAT A; 


SUPPLY  THt  NAM  OF  TH(  CURRENT  OATA  FI  LI  IN  THI  FOLLOWING 
INFILI  STATIMINT 


INFILI  '00. CUR'; 

INPUT  UStFLAG  1 

P3  TSIO  S3. 

F#  AOTCO  S3. 

•10  OCNCO  S3. 

•14  STRTOATt  MOOYYI. 
•34  OURATION 
•33  NOOSIO 
•40  NOCAO; 

FORMAT  STRTOATI  MMOOYYI . : 

IF  (OCNCO  IQ  'STO'I  THIN  00; 
OUTFUT  STOOATA; 


SUFFLY  THI  NAM  OF  THE  CURRENT  STANOARO  OATA  FILI  IN  THI 
FOLLOWING  FILI  STATIMINT 


fflWWfWf1  ■WWW1!  m  M  »  W 


FILI  ' 00. NIWHIS '  ; 

PUT  USIFLAG  1 

•3  TSIO  S3. 

9%  AOTCO  S3. 

•10  OCNCO  S3. 

•14  STRTOATI  MMOOYYI. 

«4  OURATION  7. 

•33  NOOSIO  7. 

•40  NOIAO  7.  ; 

ENO: 

ELSE  OUTFUT  TESTOATA: 

PROC  PRINT  OATA-STDOATA;  TITLE  'CURRENT  STANOARO  OATA'; 

■ 

■  COMFARt  HISTORIC  A  CURRENT  LETHALITY  RATE  (  STANOARO  OATA  I 

■  * 

OATA  COMPARE: 

MERGE  HIST  STOOATA; 

PC«N0IA0/N0Q3E0; 

RCSTAR  •  ARSIN( SQRT( PC) ) ; 

MUCJTAR  .  ARSINt SORT (MUC)); 

21»(RCSTAR-HUCSTAR)/SQRT(  3S/NOOSIO  *  SIG03 ) ; 

OSL*3»PROINORM( *1«AtS< Z1 ) ) ; 

KEEP  PC  MUC  SIG03  21  OSL; 

PROC  PRINT  NOQSS; 

TITLE  'COMPARE  CURRENT  STANOARO  TO  HISTORICAL  LETHALITY  RATE'; 
TITLES  'RATES  ARE  SIGNIFICANTLY  OIFFERENT  IF  OSL  IS  LESS  THAN  O.OS' 
OATA  CURSTO; 

MERGI  STOOATA  HIST; 

XC  •  NOE AO; 

NC  »  NOOSED; 

SIG003»SIG03; 

KEEP  XC  NC  MUC  SIG003; 

OATA  CURTEST; 

SET  TESTOATA: 

XT  »  NOE AO; 

NT  .  NOOSED; 

KEEP  XT  NT  OCNCO; 

OATA  OUMIV; 

SET  CURTEST; 

IF  _N  EG  1  THEN  SIT  CURSTO; 

SINCLUOf  CCNTAMC; 

OATA  FINAL; 

MERGE  TESTOATA  CRITICAL.  SY  OCNCO; 

PROC  PRINT  OATA«FINAL; 

TITLE  'CURRENT  TEST  OATA.  CRITICAL  NUMER  ANO  OSL'- 
TITLE2  'TEST  OECON  FAILS  IF  NOEAO  IS  GREATER  THAN  CRITNUM' • 
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DATA  I  TIN  TIMNAANT; 

srr  outwv; 

•  INCLUOS  NOMINAL  VALUIS 

'll**  •  0.10; 

SXSOO  •  SQNT( SI8002  ) ; 

SXQ01  •  0.4; 

• 

>  COMAUTI  CONSTANTS 
*  ! 

AX  »  3*AMSXN( 1.0): 

MUCSTAN  •  AMSXN(  SCAT ( MUC ) ) ; 

SI 0C3  •  . 2S/NC; 

SX0T2  •  .2S/NT; 

SXQOOa  •  SXQOO*«3; 

sxeoia  •  sxooi**a: 

wo  ■  sxoooa  /  < sxoooa  *  sxacai; 

wi  »  sxooia  /  (Sxooia  ♦  sxocai; 

AC  •  XC/24 ; 

AC ST AM  •  AASXNf  JQAT ( AC ) ) ; 

A  •  tXA(-(ACSTAA-MUCSTAA)*.2/(2«(SI0C2+Sla012))  ) 
/S0AT(2«AI«(Siac2+SXM12) )  / 
(IXM(-(ACSTAM-»tJCSTAM)»«a/(2«(SISCa*S:Q0.a) )  ) 
/S0AT(2*AI.(SI0Ca*JI0002) )) ; 

IASTAM  •  IAS»N/( 1-|AS*IAS»A) ; 

J  ' 

•  COMAUTI  OSL  AON  OSSKAVIO  TIST  DATA,  XT 
' ATSTAA  •  AASIN< SQNTl XT/NT )) ; 

ZO  •  (ATSTAA  -  (MUCSTAM*WQ«(ACSTAA-MUCSTAN)n/SQMT(SXaT2*W0>SISC2); 
Z1  •  (ATSTAA  -  ( MUCSTAA*W 1 ■ I ACSTAA-MUCSTAN ) ) ) / SQAT( SXQT2.W 1 *SZGC2 ) ; 
AXQSL  •  ( 1 -IASTAM)  •  AAOSNOAM(ZO)  *  IASTAM  «  AAOSNOAMf Z 1 ) ; 

• 

•  COMAUTI  INITIAL  QUISS 

• 

KNIW  •  MUCSTAN  ♦  (  1 -IASTAM). WO* (ACSTAN-MUCST AN) 

*  IASTAA«W1«( ACSTAM-NUCSTAM )  • 

1.64S*<  ( 1 -IASTAM ).SOAT(SiaTa*WO«SIOCa)  * 

!ASTAM«SOMT( SXST2*W1*SXGC2 )  ); 

I A  N  10  1  THIN  OUTAUT  TIMA A ANT 
OUTHI?  ITIA: 

PAOC  AAINT  OATA.TIMAAANT  NOOiS ; 

TITLI  ' AAMAMITIAS  AON  CQNT AMINATEO  NUAMAL  DXSTAXIUTION' ; 

VAN  IAS  MUC  SIOOO  SI001  WO  WI; 

XXNCLUOI  CCNTAMS: 

OATIONS  NOSOUACta; 

•4INCLU0I  CCNTAMS ; 

XZNCUlOf  CCNTAMS: 

XZNCLUOI  CCNTAMS; 

OATA  CAITICAL; 

SIT  ITIA; 

I A ( KNIW  LT  0)  THIN  KNIW«0; 

ILSI  l A( KNIW  OT  MI/2 )  THIN  KNIW  «  AI/3; 

CAITK  «  SIN(KNIW)««3; 

CAITNUM  •  CIIL1 NT  •  CAITK); 

OSL" 1.0  -  AXOSL; 

PUT  51000*  SI001.  IAS"  NC«  NT*  XC«  CAITNUM*  OSL*; 

KlIA  OCNCO  CAITNUM  OSL; 
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B-5  PROGRAM  LISTING  FOR  "CONTAMS" 
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OAT*  ITIA; 

SIT  ITIA; 

KOLO  ■  KNIW: 

ZO  •  l KOLO  -  <  »MCSTAA*WO» ( HCSTAH-MUCSTAH ) ) I  /SQAT(  SI<JT2*WO»SXOC3 1 ; 

Z1  •  (KOLO  -  (  NUCSTAA*W1«<  ACSTAA-NUCSTAAI  )  )/SQ*T(  SI0T2*W1 «SXOC2  )  : 

AX  •  (l-(ASTAA)  •  AAOiNOWMZO)  •  (AST AN  ■  AAOiNOAMl  Z1 >  *  O.St; 

FAX  •  < l-(ASTAA)  • 

(XA< -ZO«»2/2)/SQ«T<a«AI«<SXaT2*WO«SXOC2) >  * 

(ASIA*  • 

exA< -zi«-2/a)/soAT(a«Ax«(s:oTj»t<i»s:oc2)): 

KNEW  ■  KOLO  -  AX/ FAX; 

POXFF  •  100»(  KNIW -KOLO ) /KOLO; 

AUT  'XC  •  '  XC  'KNIW  •  '  KNIW  '  KOLO  •  '  KOLO  '  AOXAA  •  '  AOXAA; 
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B-6  PROGRAM  LISTING  FOR  "CRITX 
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DATA  tret: 

XNPXLE  'CRITX. OAT'; 

INPUT  NC  NT  MUC  EPS  SIQOO  SI001 : 

PILE  'STSSOUTPUT'; 

PUT  _ALL_* ; 

• 

>  COMPUTE  CONSTANTS 
•  ; 

PI  •  2*ARSXN( 1.0): 

MUCXTAR  •  AASXN<SORT!MUC)>: 

.  SXQCa  «  . 2S/NC; 

SI0T2  •  . 2S/NT ; 

SI 000 2  *  SIQOO" *2: 

SXQ012  "  SIC01»»2; 

WO  •  SI0002  /  ( SI 0002  ♦  SX0C2 ) ; 

W1  ■  SXQ012  /  ( 5X0012  ♦  SX8C2); 

00  XC  •  0  TO  24; 

PC  •  XC/NC: 

RCSTAR  •  ARSXNt  SORT! PC) ) ; 

R  •  (XP( *( RCSTAR'MJCSTAR ) «"2/( 2*( SIGC2+SIQ012 ) )  ) 
/SQPT( 2>PI*( SI0C2+SI0012 ) )  / 

( IXP( -( RCSTAR -MUCSTAR ) "*2/ ( 2* ( SIGC2+SIQ002 )  )  ) 
/SORT ( 2»PI"(  SISC2*SIQ002  >  ) ) ; 

EPSTAR  •  IPS"R/ ( 1 -EPS"£PS*R ) : 

■ 

•  COMPUTE  INITIAL  SUESS 

• 

’knew  •  MUCSTAR  *  ( 1 -EPSTAR) "WO* (RCSTAR-WJCSTAR) 

♦  EPSTAR"W1"(PCSTAR-MUCSTAR)  ♦ 

1 . 84S"(  ( 1 -EPSTAR) »SQRT( SI8T2"WO*SIGC2 )  ♦ 
EPSTAR "SORT ( SIGT2+W1 "SI0C2  )  ); 

OUTPUT; 

END; 

XINCLUOE  CONTAMS; 

OPTIONS  N0S0URCE2; 

XINCLUOE  CONTAMS; 

XINCLUOE  CONTAMS : 

XINCLUOE  CONTAMS; 

OATA; 

SET  ITER; 

IP (KNEW  LT  0)  THEN  KNEW"0; 

ELSE  I P ( KNEW  OT  PI/2)  THEN  KNEW  »  PI/2; 

CPITK  •  SIN( KNEW )"*2 ; 

CPXTX  •  CEILf NT  «  CPITK); 

PILE  'SYSSOUTPUT'; 

PUT  XC*  CRITX*  ; 

PPOC  PRINT; 

VAR  NC  NT  MUC  EPS  SIGOO  SIQ01  XC  CRITX; 

TITLE  'INPUT  PARAMETERS  AMO  CRITICAL  VALUES  POP  XT' 


